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ABSTRACT 




The problem of imaging a planetary surface is considered from the 
point of view of continuous wave bistatic-radar in which transmissions 
originating on the Earth (or on a spacecraft) would be received on a 
spacecraft (or on the Earth) after reflection from a planetary surface. 

An electromagnetic wave reflected by a planet may be considered as 
a superposition of the waves from elementary scatterers on the surface. 
Using a superposition model the ordinary radar brightness distribution, 
the differential radar cross-section, the polar scattering diagram, and 
the radar albedo of a surface are all related measures of the local sur- 
face properties. 

The surface brightness distribution may be determined by cross- 
correlating the scattered fields (as measured by a spacecraft along some 
fraction of its trajectory) with the expected signal from each point on 
the surface. The predicted azimuthal resolution in wavelengths is in- 
versely proportional to the angle subtended at the target point by the 
fraction of the trajectory over which the data are taken; in range the 
resolution in wavelengths is inversely proportional to the square of the 
same angle. The feasibility of the method depends on the use of the il- 
luminating wave on board the spacecraft as a frequency reference to achieve 
the requisite stability. 

Physical analogs of the process exist as modifications of holograms 
and synthetic antenna arrays. An additional analog is that of a bank of 
tracking filters, each of which is adjusted to receive the signals from 
a separate portion of the surface. 

The maximum-likelihood estimator for the brightness of a specific 
scattering area, in the presence of white Gaussian noise at low input 
signal -to -noise ratios, is a Hilbert quadratic form with the expected 
autocorrelation function of the signal from the scattering area as a 
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kernel. For maximum resolution (minimum area) this is equivalent to the 
cross-correlation method for obtaining the brightness distribution. Fur- 
thermore, this estimate is efficient in the sense of being unbiased and 
of having minimum variance, for the low signal-to-noise ratio case. 

The estimator may be realized with well-known forms of time-varying 
or time-invariant filters, or with correlators. A new realization is 
that of a hologram scanned with a properly weighted illuminating wave. 

It is concluded that bistatic-radars operating between a ground 
station and a vehicle in space constitute a powerful technique with which 
to map and to study planetary surfaces. 
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I • INTRODUCTION 
A. STIMULUS FOR THIS RESEARCH 

The first radars, dating from experiments at the turn of the century 

1 2 * 

were bistatic in their operation. ' That is, they were radars in which 

the transmitter and receiver occupied not the same, but two separate, 

remote, and generally fixed positions. Energy was continuously radiated 

from the transmitting site, and target detection relied upon scattering 

enough of the transmitted signal power into the receiving antenna by the 

target to overcome the noise of the receiving system. The receiving 

sites were usually effectively shielded from the transmitter, and target 

detection was signaled by an increase in the power output of the receiver. 

An alternative was to allow a small amount of power from the transmitter 

to leak directly into the receiver. In this second case a Doppler- 

shifted signal from a moving target appeared as a beat note at the Doppler 

frequency. By the early 1930s, enough progress had been made to permit 

the detection of aircraft at ranges of a few tens of miles. 1 

Now, three decades later, astronomical bistatic radars, employing 

transmissions between a station on the ground and a receiver on board a 

spacecraft in interplanetary space, are being used as a scientific tool 

to probe the electromagnetic properties of the intervening medium. These 

experiments are characterized by the continuous presence at the receiver 

of a strong signal, directly from the transmitter, which is examined for 

perturbations in phase, amplitude, or polarization introduced by material 

3 

along the path followed by the wave. Independent measurements of two 
or more of these signal characteristics made simultaneously on several 
frequencies may be employed to magnify the effects of the medium and to 

3 

calibrate the experiment. 

The concept of bistatic-radar for probing tenuous astronomical media 

was quickly extended to include proposed experiments for determining the 

4 5 

properties of planetary atmospheres and ionospheres. * Here the experi- 
menter makes use of his control over the spacecraft trajectory. By causing 
the spacecraft to pass behind the planet, the propagation path between the 

•if 

References are listed at the end of the report. 
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ground station and the vehicle is forced to slice through whatever sur- 
rounding medium there may be. Since the spacecraft's position with 
respect to the planet is well known, it is possible to associate the mea- 
surements with particular paths through the medium and consequently with 
particular portions of the atmosphere. The data obtained, together with 
certain auxiliary information, are then used to establish the salient 

physical and chemical properties of the atmosphere and ionosphere in- 
6,7 

volved . 

The next possibility is immediately clear: A spacecraft may serve 

as one terminal of a bistatic-radar system for exploring the surface of 
a planet. Either transmissions from the earth would be directed at the 
planet and the scattered energy would be received on board the vehicle, 
or energy directed toward the planet by the spacecraft would be received 
on the ground after scattering from the surface. There are several rea- 
sons to believe this will be useful. 

1. Relative to purely passive observations, the probing signal is 
under the experimenter's control, and its characteristics may 
be chosen to maximize the experimental sensitivity to a given 
surface property. 

2. Relative to ground-based radars, (a) one leg of the radar path 
may be made very short, thus increasing the signal strength at 
the receiver, and (b) the angles of incidence and reflection 
are variable, thus allowing the resolution of the surface-scat- 
tering properties as a function of direction. 

3. Relative to on-board radars, only a transmitter or a receiver 

is required thus considerably reducing the weight, power con- 
sumption, and complexity of the spaceborne equipment. If this 
advantage is to be realized, then only one-way transmissions may 
be used for the experiment: signals may go from the spacecraft 

to the ground or vice versa. Data stored on the spacecraft will 
be retrievable by telemetry. 

4. Relative to optical measurements, radar may be the only means 
for observing the surface of Venus for some time. 

5. Signals reflected from the surface may be needed to remotely 
probe the atmospheres of some planets below a critical refraction 
level . 


SEL-67-042 


2 


*4 


This research was instigated to determine if there were ways in 
which the bistatic technique could be used, extended, or modified to 
realize the potentials listed above for the study of planetary surfaces. 

B. PREVIOUS CONTRIBUTIONS 

It is not clear who first suggested the bistatic-radar approach to 
planetary surface studies. To this author's knowledge, the first written 

G 

comments were in a January 1961 memorandum, by Evans and Pettengill, in 
which they discussed the limitations of monostatic radar for surface 
studies and means of resolving these difficulties with bistatic-radar 
probes. Their suggestion was to transmit with an omni-directional antenna 
from a vehicle orbiting the Moon and to receive the scattered energy with 
a tunable radiometer. In effect, they would rely upon strong quasi-spec- 
ular scattering to isolate a fraction of the disk and resolve the mono- 
static scattering law ambiguities. 

The second written contribution evidently came from a NASA study 
group convened for the purpose of discussing bistatic-radar, with the 

result that the potential advantages listed in the introduction were 

9 

articulated in the group's report, which was authored by Eshleman. What 

is important here is that the application of bistatic-radar was realized 

to be a versatile and potentially powerful tool for planetary exploration. 

Finally, in a work dealing primarily with bistatic probing of plane- 

4 

tary atmospheres, Fjeldbo derived the properties of the quasi-specular 
scattering from a surface as seen from a space probe, and related the 
spectral characteristics of the scattered signal to the parameters of a 
surface generated by a stationary Gaussian process. 

There have also been several papers and reports on the degradation 
of communications systems caused by reflections from planetary surfaces 
by extrapolation of monostatic radar results. 

In addition to the work above, there were a number of meetings and 
conferences 1 ^, beginning in June 1960, when the general topic of bistatic- 
radar for astronomical studies was introduced by V. R. Eshleman. Several 
proposals for studies of the use of lunar orbiters and planetary flybys 
were initiated at Stanford Univeristy in 1962. 
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The first bistatic-radar experiments for the study of the surface 
of a celestial body were carried out in October 1966, by the author, in 
collaboration with his colleagues at the Center for Radar Astronomy, 
Stanford University . ^ 

C. DESCRIPTION OF PRESENT WORK 

The present work proceeds along the following lines. 

Chapter II is a general examination of the scattering properties of 
a planet as viewed from a spacecraft. The analytic signal representation 
of the scattered wave is introduced, a simple integral representation 
of the scattered wave is derived in terms of the analytic signal, and the 
radar brightness distribution is defined. The connection between this 
formulation and the more common radar cross-section and polar scattering 
diagram is given next. Specific features of the scattered fields as 
received on board a spacecraft are given last. 

Chapter III deals with a method for recovering the radar brightness 
distribution from the scattered wave on the basis of a heuristic argument. 
The resolution cell, signal-to-noise ratio, and sensitivity to errors are 
determined for this method. It is also shown how the method may be made 
self -ca libra ting by use of the illuminating wave as a reference signal. 

Chapter IV gives three physical interpretations of the inversion 
method. First, it is considered as a synthetic array, performing as an 
antenna focused in its own near-field. Second, we give an interpretation 
as a tracking filter and explain the resolution properties in terms of 
the filter characteristics. Third, we develop an analogy with holograms, 
and show how the inversion may be automatically realized with a properly 
constructed diffracting screen. 

Chapter V develops the inversion as a maximum-likelihood estimator 
of the brightness distribution from a specified surface area. The general 
form of the estimator is derived. Then after restricting the problem to 
low input signal-to-noise ratios, the estimator is applied to the inver- 
sion problem. It is shown that the process produces an unbiased minimum 
variance estimate of the surface brightness distribution. The relation- 
ship between the variance of the estimate of the surface brightness and 
the size and shape of our resolution cell is investigated in the last 
section . 
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Chapter VI presents five methods of realizing the estimator. Realiza- 
tion in terms of two time-invariant filters (one physically realizable 
and one not physically realizable), a time-varying filter, and an autocorre- 
lator are given. The realization of the optimum filter through the holo- 
gram analog is also discussed. 

Chapter VII gives some experimental results from the first applica- 
tion of bistatic-radar to the surface of a celestial body. An experiment 
using Lunar Orbiter spacecraft is described and some samples of the obser- 
vations are given. A two-dimensional map, made from observation of a 
reflected, continuous-wave signal, is presented. 

Chapter VIII gives the conclusions we might reach from the preceding 
seven chapters. 

The author has made the following contributions to this research: 

1. A method for the two-dimensional bistatic-radar mapping of 
planetary surfaces using a continuous -wave mode of transmission 
between the Earth and a spacecraft. 

2. An analysis of that method based upon techniques from conventional 
radar and optical imaging theory. 

3. Three analog interpretations of the method in terms familiar 
devices . 

4. The derivation of a specific maximum-likelihood estimator for 
the radar brightness of a general scattering area from data 
taken on a moving spacecraft. The conditions under which the 
original method is optimum and the means by which it should be 
modified for other conditions are determined. 

5. The first bistatic-radar measurements of the scattering proper- 
ties of the Moon, and the preliminary reduction of the data in 
terms of the ideas developed above . 
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II . PHYSICAL ASSUMPTIONS AND DESCRIPTION OF THE SCATTERED WAVE 
A. THE ANALYTIC SIGNAL 

The use of complex numbers to describe real physical processes is 
a familiar tool of scientists and engineers. Sinusoidal analysis in the 
time domain frequently makes use of the relations 

cos fit = R^{e'^ t } 

and (2.l) 

sin fit = I {e^^} 
nr 

to interpret a real oscillation as the projection of a rotating phasor 
onto a fixed axis. Similarly, complex Fourier analysis exploits 

jfit -jfit 
e + e 

cos fit = 


and (2.2) 

jfit -jfit 
. ^ " © 

sm fit = 

^ J 

to decompose real functions into complex integrals over frequency. 

While these formulations have great merit, both in a mathematical 
sense and in the physical insight they afford, they are somewhat awkward 
for dealing with signals when time-frequency dependence is important. 

In this report we will be dealing extensively with time changes in 
relative phase among a group of signals. Consequently, the analytic signal, 
a formalized, time-varying version of the rotating phasor, is introduced. 
This device will allow us to deal explicitly with the relative phases of 
signals, and the time variations of phases. The presentation here is 
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standard, although in what follows we shall try to emphasize the fea- 
tures of this representation which will be of use to us later. A few 
elementary results from communication theory will be stated as they apply. 

We shall consider a real signal s(t) with Fourier spectrum s(o>), 
so that 

r 00 

s(cd) = J s(t) e dt 


and 


(2.3) 



last 

e 


do> 


As a direct consequence of s(t) real we have 

S(no)) = S(o))* (2.4) 

where the asterisk is used to denote the complex conjugate. Hence, nothing 
is lost if we consider only s(cjd) for a) > 0, since the function of the 
negative argument may always be recovered from its value along the posi- 
tive real axis. 

The time function generated by transforming the truncated spectrum 
obtained by discarding s(co) for 00 < 0 and doubling the remainder is 
the analytic signal 


z(t 




do) = x(t) + jy(t) 


(2.5) 


The function z(t) is not real, but its real part, x(t), is the original 
s(t ) , as may be shown by 
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dcu = 


s(t) 



Jcot 

e 


*(*) + **(<) , R f (z(t)) (2.7) 


or s(t) = x(t). 

Eq. 2.5 allows us to interpret z(t) as a collection of infinitesimal 
phasors, the e^^dco, rotating in the same direction, albeit at different 
OD, with respective amplitudes s(oo). If s(t) is considered as a narrow- 
band signal with spectrum s(cd) which is non-zero only about some fre- 
quency £7, this interpretation becomes quite clear. For example, take 


s(a)) = f (o)-ft) + f (“tjd-£7), f(cn) = 0 for |cjd | > W 

W < ft (2.8) 


as a representation of a narrow-band signal. The analytic signal becomes 


z(t) 



j (co+n)t 

e 


dcu 


jilt 1 
e — 

Tt 


f f(a>) dcu = Z(t) 

-W 


(2.9) 


As before, 


s ( t ) = R i {z(t)) = R^{z(t) e jat ) (2.10) 

Writing z(t) = x(t) + j Y(t), we have 

s ( t ) = X (t ) cos ftt - Y( t ) sin fit (2.1l) 

The x(t) and Y(t) are referred to as the in-phase and quadrature com- 
ponents of s(t). The complex envelope of s(t) is z(t). It should be 
noted that ft has been specified in only the most general terms, and that 
the form of these results is not critically dependent on the value of ft, 
even though the exact z(t) obtained may be. In addition, z(t), and 
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X(t) and Y(t) likewise, are relatively slowly varying functions that 
change significantly only in a time greater than l/w. 

In terms of Eq. 2.5 we have separated z(t) into a reference phasor 
or carrier of constant frequency and amplitude, and a complex 

modulating function z(t) of amplitude 

m( t ) = {X 2 (t) + Y 2 (t)} /2 


and angle 


pit) = tan" 1 (2.12) 

From now on the term phase will refer specifically to a function pit); 
amplitude will refer to M(t). In particular if M(t) is constant, so 
that we have a purely phase-modulated wave, we picture the signal as a 
phasor spinning in the complex plane with the nominal angular position 
(phase) fit and gaining or losing by the quantity pit). The represen- 
tation of cos ( fit - 0(t)) is simply e^ fit ~ Thls con cept will 

be important to us later. 

A distinct virtue of this formulation is the ease with which x(t) 
and Y( t ) may be obtained by measurement. If we multiply s(t) by 
cos fit the result is 

s ( t ) cos fit = i x(t) + i x(t) cos 2QX - ^ Y(t) sin 2ftt (2.13) 

Removing the 2flt terms by filtering and doubling the resultant gives 
X(t). A similar procedure employing sin fit yields Y(t), A block dia- 
gram for this operation is given in Fig. 1. The cutoff frequency for the 
filter, o> c , need only be low enough to avoid the twice-f requency terms 
in Eq. 2.13, and high enough to pass all of X(t) and Y(t). Generally, 
it is advantageous to use the lowest possible value consistent with the 
modulation bandwidth. 
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FIG. 1. METHOD FOR OBTAINING COMPLEX ENVELOPE. 

We may also relate z(t) to the energy E in the signal. From 
Parseval 1 s theorem 


f +°° 0 /*+°° 9 

|s(oj)| da> = J s (t) dt ( 2 ‘ 14 ) 


f |s((n)| 2 d6D=2 f | s (c 

■—*00 * / 0 


<&d for S 


(-0)) = S*(co) (2.15; 


and since x(t) = s(t) we have 


r -k>o _ /* 4 oo r+°° 

= f x 2 (t) dt = f Y 2 (t) dt = - J |z( 

-Lv, •'-00 ^-OO 


t )| 2 dt (2.I6 


If we now consider filters with real impulse response h(T), and 
transfer functions h(oj), it follows that we may define a complex im- 


pulse response 


e(t) =~ f°° h(co) e^ T dco 

^ Jc\ 


(2.17 
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As before, we have generated the complex representations by discarding 
the spectrum for co < 0 and doubling the resultant. By restricting the 
problems to narrow-band filters about ft, we again have 

g(x) = g(t) e jfiT (2.18) 

with G(t) the complex envelope function of the filter response in a 
manner completely analogous to z(t) and Z(t). Still with the narrow- 
band restriction and letting s^t) and s Q (t) be the input and output 
signals respectively of h(t). 


/ 00 

h(T) s (t-T) dT 

no 


(2.19) 


It is a simple matter to show that 


/•+00 

Z 0 (t) = e jat I G(t) Z (t-T) dT (2.20) 

•'-oo 

where Z (t) is ihc complex envelope of the analytic response function 
of the filter. Thus, the complex envelope of the response of a narrow- 
band filter to a narrow-band signal may be calculated directly from the 
complex envelopes of the input signal and the impulse response function 
of the filter. It is also possible to apply the complex description of 
signals to band-limited noise waveforms. We define N(t) as the complex 
envelope of a narrow-band noise waveform n(t). Thus 

n(t) = Rg{N(t) e^} 

as before and 


N(t) = X (t) + j Y (t) 
n n 


( 2 . 21 ) 
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defined the in-phase and quadrature components of the noise. Clearly, 
these two components of the noise may be obtained by sine and cosine 
multiplication of n(t), followed by filtering, exactly as before. 

Letting $(a>) represent the power spectrum associated with n(t), 
the autocorrelation function of the process is 



( 2 . 22 ) 


The corresponding analytic function c(t) is given in this case by 

J r 00 

®(o>) da; = c(t) e JflT (2.23) 

0 

13 

A well-known result of this is that for stationary narrow-band processes 

< X „ (t l ) X „(‘2)> = < Y Jh ] V „ (t 2 ) >= 

(2.24) 

< Y „ (t l ) X „ (t 2>> = < X „ (t l ) Y „ (t 2>> = I m< C(t l- t 2 ) > 

C(t -t ) is the complex envelope of C(t). We shall make use of two 
-L £ 

consequences of this result: 

1. If the noise is a Gaussian process, the X (t) and Y (t) 

n n 

are also Gaussian since they are the results of linear operations 
on Gaussian random variables. 

2. Since 0(a)) must be real, c(t) must be of the form 

Job T 

C (t) = {e(t) + j 0 (t) ) e c (2.25) 

where e(t) and o(x) represent even and odd functions respectively. 
If, in addition, 0 (cn * t ft ) is even about the co 1 origin (perhaps 
requiring a suitable choice of ft), then in addition o(t) = 0. 
c(x) = e(t) + j(o) and thus, x n ^ t ) and are stat istieally 

independent from Eqs. 2.24. 
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B. THE SCATTERING PROCESS 


If we now consider the fields in a common reference frame originating 
from an oscillator at frequency ft, the Maxwell equations in homogeneous 
media are 


V X E = -jftpH 

V x H = jfteE + J 

(2.26) 

V • H = 0 

V • E = p/e 

We have used our complex representation to write vectors as V e J ^ t . 

Any divergenceless vector is the curl of some other vector, so 

H = V X A (2.27) 

Substituting the above into Eqs. 2.26, we have 

V X (e + jftA) = 0 (2.28) 

where A is the magnetic vector potential. 

Now, invoking the complementary relationship that a curl-free vector 
is the gradient of some scalar gives 


or 


E + jftA = -V0 
E = -V0 - jftA 


(2.29) 


The function 0 is the electric scalar potential. 

Choosing V • A = - j fie 0 , so that A corresponds to a Lorentz 
gauge, 
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and 


(2.30) 


E = -JAA + V(V • A) 


H = V X A 

We shall use the expression for A in terms of the volume current 
distribution J (r ) . 


ig) -a /, 


J(r') 


— jk | r-r ' 


dv' 


(2.31) 


-t*' 


r-r 


Here, r and r 1 are position vectors used to denote field points and 
source points, respectively, and to emphasize that the A, to be evalu- 
ated at r, results from j(r') integrated over V 1 . As usual k is 
the magnitude of the wave vector 




A general illustration is given in Fig. 2. 



FIG. 2. GEOMETRY FOR EQ. 2.31. 
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Reintroducing the e 


jftt 

» 


A (?) e Jftt = JL f J(r')e~ jk i r - r ' 

4jt ./ , 

* i j * 




dv 1 


( 2 . 32 ) 


r-r 


we see that A(r) is the complex envelope associated with out representa- 
tion of the magnetic vector potential. Although at the moment A is 
a constant in time, it will not remain so. The electric field will be 
determined from Eq. 2.30 by substituting the representation of 2.31 
for the vector potential. 

Calculating V- . A(r) 


V • A (r ) = -jp- V- 
r ' 4 rt r 


L 


J(r') 


-jkl r-r' 


dv 1 


— T* * 


r-r 


(2.33) 


= zJilii 

4jt 



rr 1 


dv 1 




u — 
rr 1 


-jk Ir-r'l 




dv 1 


r-r 


(2.34) 


u — 
rr 


1 


— T* 1 


r-r 


— T* * 


r-r 


By restricting ourselves to r far removed from V f , the second term 
may be neglected with respect to the first. Taking V-(V- • A) gives 


v-(v- • a) = - ^ r 

r r 4jt J , 


(J(?«) • u~ f )u~, 
' rr f rr 1 


dv 1 




r-r 




l 


(j(r') 


u . )u — e 
rr' y rr 


-jkl r-r 1 




dv 1 


(2.35) 
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Dropping the term in l/|r-r'| again and substituting in Eq. 2,30 

-jkl r-r ' 


•Cr) - -3 % 

v r-r 1 


-jk r-r 1 


(j(r') -u )u e 


. Qu f ' ' ' rr ' rr ■ 

av + i, ~ 


dv 


r-r ' 




= -j 


4n 


/, 

i j » 


(j(r') ' u--, )u - J(r' ) e 

rr rr ' 


-jk |r-r f | 


dv' (2.36) 


r-r ' 


The bracketed terms in the integral are the components of j(r*) normal 
to the direction of propagation. This will be denoted by 


J (r 1 ) = (j(r f ) ’ u — ( ) u — f “ J(r f ) (2.37) 

n rr 1 rr f 

By ignoring the higher order terms in the demoninator and assuming r 
far removed from the volume V 1 we have restricted ourselves to the 
radiation fields. Except when explicitly noted this will be the case 
from now on. 

Equation 2.36 illustrates an important property of electromagnetic 
radiation; i.e., except for a constant, the radiated fields in any direc- 
tion may be considered as a linear superposition of the currents that 
produced them, retarded in phase by k|r-r f |, and attenuated by l/|r-r*| 
In particular, E(r) results from the convolution of a point radiation 
source, e’ jk l r l/|r| and the function J^r* )• Furthermore, the only 
physically restrictive assumption about the nature of 
it exists in electrically homogeneous and unbound media. 

These results apply almost automatically to the radiation from plane- 
tary surfaces. The volume over which our currents exist becomes a thin 
region containing the part of the planet in which we are interested, 
while the currents themselves tend to sheets. For the surface of a planet 
we might have a situation similar to that given in Fig. 3. Thus we write 

J r £7-f\ -jklr-r'l 

f Bill e ! L ds » (2.38) 

c f I ^-r 1 I 
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when S' = the surface of integration and K(r') = current sheet on s: 
We do this with the understanding that K(r* ) may in fact have depth, 
and that if it does, we must take it into account by reverting to inte- 
grals over volumes. 

In more complicated situations, where the fields can no longer be 

considered monochromatic, it is possible to extend the formulation that 

14 15 

we have outlined directly via the Fourier transform. * However, that 
will not be necessary here. 

To use these relations in the problem at hand, we adopt the view 
that the currents in our integrals are the result of incident radiation 
from a distant source impinging upon our surface. The reradiated, or 
scattered waves, are then given by Eqs. 2.37, 2.38. Unfortunately, the 
relationship between these currents and the source fields producing them 
is one of the most difficult problems in electromagnetic theory, and one 
in which little progress has been made. Fortunately for us, however, 
the form of the scattered wave given is sufficient for our work. The 
important property of the scattering, and the one that is critical to 
the later developments in this study, is that no matter what the surface 
is like, or how the surface currents are generated, the radiation fields 
can always be considered as arising from a continuous distribution of 
point sources. 


C. THE POLAR SCATTERING DIAGRAM AND BRIGHTNESS DISTRIBUTION 

Most radar measurements ultimately consist of a determination of 
the power scattered into the receiver by the target for some particular 
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set of antenna and receiver parameters. The quantities of interest might 
be the signal level present on an absolute scale, the ratio of the trans- 
mitted-to-received signal, or both the absolute and relative values. Con- 
sequently, it is necessary to have a parameter that describes the target 
in terms of a power measurement by the radar. Such a parameter obviously 
serves two purposes. The first of these is to provide a standard for 
comparing different radar systems. The second, which is particularly 
applicable to radar as a scientific probe, is to characterize the target 
itself. The quantity apropos this usage is the target radar cross-section, 
which is defined as the capture area required to intercept exactly that 
amount of power from the transmitter, so that if it were to be isotropi- 
cally reradiated, it would scatter the same power into the receiver as 
does the actual target. Using this definition, the radar equation is 
commonly written as 


P R = 


PG 
T T 

2 

4jtR, 


4*R, 


R 


(2.39) 


where P and P are the transmitted and received powers; R and R 
T R 14 

are the ranges from the transmitter to the target and the target to the 
receiver, respectively; A is the effective area of the receiving an- 
tenna viewed from the target while is the gain of the transmitting 

antenna in the target direction; cr is the radar cross-section that we 
have just defined. A sketch showing these quantities is given in 
Fig. 4. 

When the target is very large compared with the other dimensions 
involved as in Fig. 5, the usage must be modified slightly. Since in 
this case it is usually possible to distinguish one part of the target 
surface from another, a differential radar cross-section is introduced. 
Thus , we have 



(2.40) 
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FIG. 4. ILLUSTRATING THE RADAR EQUATION. 

where <7 is the number of square meters of radar cross-section per 
unit area of target surface. The radar equation becomes 


% 


PGA <7 
T T R 

/ .222 

(4rt) RjRg 


As 


(2.41) 



The term AP is the power scattered by the differential element As. 
o R 

a has immediate expression as a ratio of the scattered and incident 
fields introduced in the previous section: 
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(2.42) 


It is clear from even the most cursory examination of Eqs . 2.37 and 
2.38 that cr° is, in general, a function of both the direction of inci- 
dence and the direction of reflection. Introducing the polar angles of 
Fig. 6 , 0 ^, 0 ^, 0 as the angle of incidence, the angle of reflection, 

and the angle between the planes of incidence and the plane of reflection, 
respectively, allows explicit notation 


= o°(e . e ; 9.) 


(2.43) 


The surface itself is assumed homogeneous in its scattering properties, 

so that a° is independent of rotation about the normal from which 9^ 

and 0 are measured, 
r 


- „ \ 








FIG. 6 . SCATTERING GEOMETRY. 


The polar scattering diagram ^(0 r > 9p> 9^) norma ii zec * radar 


cross-section 


^ 9 r’ 9 p : 9 i ) " /* 2 « /-jt /2 


a °(0 , e ; e ± ) 


I I a°(e r , 0 ; e.) sin e r de r d0 r 

* / o J o 


(2.44) 
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shown in a manner completely analogous to the directivity of an antenna. 
Letting the total reflectivity be 

X 2t t /*rt/2 

J <r°(e p » e p ; 9.) sin e r de r de p (2.45) 

we write 

cr° = p(e ± ) *(e r , 9 p ; 0 i ) (2.46) 

Thus, a determination of o°(0 r , 0 p ! 9^ is equivalent to a determina- 
tion of the polar scattering diagram and a reflectivity factor p(0,). 

While radar cross-section and the polar scattering diagram are 
extremely useful for describing the scattering from a single target or 
from an element of surface, they become awkward to use when expanded to 
include surfaces that exhibit considerable variation in their scattering 
properties from one point to another. 

To escape these difficulties, we adopt another point of view. If 
we suppose for the moment that the surface is examined from some fixed 
vantage point, so that our angles of incidence and reflection are fixed 
for each portion of the surface, then we may introduce a brightness 
function B(r), which depends only upon surface location. Using the 
standard definitions of total brightness 1 ^ and cr° as the differential 
radar cross-section at r for our particular angles of incidence and 
reflection, we have 


B(r; 0 > 0 » 9. ) 
r p i 



4jtR 


(2.47) 


Assuming that the incident power per unit area of incoming wave is 
lw/m^ gives 


PG 
R T 


4*R 


1 
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and 


B(r; e , 9 , 0 ) 

r p i 


a °^ ; 9 r ’ V 0 i ) 

4 it 


(2.48) 


Thus, in effect, we have three names for the same thing. The brightness 
distribution B(r), which we use when describing the power arriving 
from various portions of our surface; the polar scattering diagram 
which we use when describing the angular dependence of the scatter of 
an elementary surface element; and the radar cross-section <j°, which 
we use for convenience in the radar equation. Each of these has hidden 
in it the parameters and dependencies displayed by the other. 

The albedo of a surface is the total brightness when integrated 
over the hemisphere above the surface for a given angle of incidence. 
From Eqs. 2.48 and Eq. 2.45, this is related to the reflectivity by 


Albedo = 



(2.49) 


Assuming K(r) is the result of E i , the brightness may be easily 
related to the current distribution on the surface. Differentiating 
Eq. 2.38 with respect to s' we find 


AE 


4 * 


K (r 1 
n 


-jkl r-r 1 


i r-r ’ 


(2.50) 


so that 


and 




2 


B(r') a |K n (r')| 2 


(2.51) 
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The brightness distribution is proportional to the magnitude squared of 
the current distribution on the surface. 

D. THE SCATTERED SIGNAL RECEIVED ON BOARD A SPACECRAFT 

The developments of the preceding sections can now be used to de- 
rive some rather general properties of the radiation scattered by a 
planet. Having done so, it will be possible to predict the character- 
istics of the signals which will be seen by a bistatic-radar operating 
between the Earth and a spacecraft near the planet. 

There are obviously two ways in which such a radar could operate. 
Either the transmitter remains on the ground while the receiver is 
carried into space, or the transmitter must become spaceborne while the 
receiver is left behind. Each of these two possibilities has its own 
advantages and disadvantages in practical terms. From a theoretical 
viewpoint, however, the reciprocity theorem assures us that they are 
equivalent, since the fields present at the receiver terminals must be 
the same in either case, all else being equal. But there are conceptual 
benefits to be gained from considering the first, or up-link, case where 
it is the receiver that is carried by the spacecraft; therefore, this is 
the case we will analyze. 



Consider the situation depicted in Fig. 7. The origin of the coor 

dinate system shown is at the center of the planet. In this system, 

r r 1 and r are position vectors of the Earth, a point of the sur- 
E 

face of the planet, and the spacecraft, respectively. The planet is 
illuminated with radiation from the Earth of the form 
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(2.52) 


E. = E 
1 o 

As usual k. is the incident wave vector in the direction of propaga- 
i — 

tion. The Earth is assumed to be so remote that the magnitude of E^ 

is constant near the origin of our coordinate system. The resulting 

scattered wave is denoted E which, with the representation of Eq. 2.38, 

s 

may be written as 


jCftt-k^rg) 

e 


e 

s 



jftt 

e 


jfit 

e E 

4* 


L 


K(r' ) 


-jkl r-r' 


1 ds 


-T- ' 


r-r 


(2.53) 


As before, we need only concern ourselves with E , the complex envelope 

~ S 

of the analytic signal representation e^. The area of integration s^ 

is the surface included within the spacecraft horizon on the planet. 

The comes directly from the assumption that the current and the 

scattered fields are caused by the e. . The field distribution E is 

1 a 

pictured as a fixed pattern attached to the planet in much the same way 
that an antenna pattern is associated with a physical antenna. 

The spacecraft antenna pattern is taken into account by introducing 
an auxiliary coordinate system, parallel to the original and centered 
on the spacecraft as in Fig. 8. In the auxiliary system, let 6 V ( U ) 
denote the voltage gain pattern of the spacecraft antenna for a wave 
arriving from the u direction. Modifying the integral by weighting 



FIG. 8. ANTENNA VOLTAGE GAIN. 


SEL-67-042 


24 



each contribution to the field by the antenna gain in the direction from 
which it arrives, the voltage at the antenna terminals except for a con- 
stant is 


V(?) = f 

J s t 


K(r') 


-jkl r-r' 


-y< 


g (u — ,) ds 

r* y*r* » 


r-r 


rr 


(2.54) 


Fjeldbo has used an expansion similar to the method of stationary 
phase to analyze the case of scattering from a planet whose surface is 
generated by a zero mean, Gaussian random process superimposed upon a 
spherical surface. Starting from the Huygens— Khirchoff approximation 
and expanding about the point of specular reflection for the mean surface, 
Fjeldbo concluded that the average radar cross-section of a gently un- 
dulating planet is equal to that of a smooth sphere of equal size, and 
that the principal portion of the energy is returned from a region about 
the instantaneous specular point. The size of this region is determined 
by the probability, as one moves away from the specular point on the sur- 
face, of a plane tangent to the local surface being properly oriented to 
produce a reflection toward the spacecraft. In effect, his method selects 
the portions of the surface that tend to produce mirror-like reflections, 
in that the phases of the current distribution on these surfaces undergo 
a linear variation with r' . Such surfaces are frequently thought of as 
plane facets or as gently curved surfaces which vary smoothly on the 
scale of a local Fresnel zone. Numerous other authors, while not applying 
their results to a spherical geometry for planets, have come to similar 
conclusions . 

Contributions to the integral from regions where the phase is not 
stationary are not so easily handled. It is generally supposed that 
these regions correspond to portions of the surface that are extremely 
rough, or inhomogeneous, or vary in some other way that, on the scale of 
a wavelength, causes large phase fluctuations in our integral Eq. 2.53. 

The reradiation associated with these rapid phase fluctuations of the 
current sheet is called the diffuse component . 
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Generally, the fraction of the field generated by the diffuse com- 
ponent is much smaller than that from the quasi-specular . However, since 
the fields associated with these two mechanisms originate on different 

parts of the planet, they have different u-- , and thus may be enhanced 

rr r 

by a suitable choice of antenna function g . If the incident wave 

v 

changes relative to the surface, as it must do when the planet rotates, 
a new pattern is generated. However, until the rotation is great enough 
to cause significant variations in the magnitude of the current distrib- 
utions, the effect is to modify the phases of the current distribution 
by the change in the phase path of the illuminating wave. This may be 
displayed in our integral by factoring out the phase retardation of the 
incoming wave. 


E 

s 



- jk • r 1 - jk I r-r ' 
K T (r T ) e 1 


-T*' 


r-r 


ds 


_ _ +jk -r’ 

K' (r' ) = K(r' ) e 


(2.55) 


Thus far, e^ has been treated as a point function observed at r. 
But recognizing that a spacecraft must be in motion and that planets 
rotate, the two position vectors r f and r must now be considered as 
functions of time. Writing this out: 


E g (r(t)) 


E_ 

4n 



-jk | r ( t ) -[r ' (t )] | ~ jk * [r ' (t ) ] 
K'([r'(t)]) e l 

|r‘(t) - [;*(t)]| 


ds 


(2.56) 


where [r(t)] is the retarded quantity r' ^t - I r ' — hj . The phase 

from any point on the surface is 


Cp(t) = -k i *[r , (t)] - ku-- f (t) * (r (t ) - [r'(t)]) (2.57) 
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The unit vector u-_,, must now be defined in terms of ?, and the 
retarded Frequency is the time derivative of phase. So 

= dt = * t r '] " k-~, • (r - [r* ] ) (2.58) 


This is simply the expression for the Doppler shift associated with the 
relative motion of the points r and r ' . The total scattered signal 
is the composite of many small signals of different frequencies. It is 
clear from Eq. 2.58 that E g (t) is band-limited, since all the terms on 
the right are finite. Furthermore, it is also easy to show that E (r(t)) 
represents a narrow-band signal for all spacecraft velocities much less 
than the speed of light. Thus, we see that the effect of spacecraft mo- 
tion is to map points on the planet into the instantaneous frequency do- 
main by the scalar product of the velocity vector and the wave vector. 

This mapping is not necessarily one-to-one, since several points may share 
the same instantaneous Doppler frequency. However, it should be clear 
that points close together on the surface must remain close in frequency. 
Bright spots on the surface will show up a strong component of the fre- 
quency spectrum of the scattered signal. 

It is also convenient to be able to describe the variation of the 
scattered fields without reference to time. For this we use path length 
along the trajectory as a coordinate with the vector relations 

s(t) = f v(t) dt (2.59) 

J t 

o 

and 

i(t) = ;(t) = jjf 

Distance measured along the path will be called £. We define the spatial 
frequency in the direction u as 


03 s V Cp . u (2.60) 
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This is easily related to the Doppler frequency by 


v 


cd = (Vcp • u)| v| = cd| V | (2.61) 

For a spacecraft less than a few planetary radii from the surface, the 
velocity of the vehicle is a great deal higher than that of a point on 
the surface, 


r' « 


so that the Doppler shift is dominated by the spacecraft velocity 


cd = k — . • u- 

rr 1 v 


where 

u- = ;/|v| 

With this approximation, the spatial frequency is just the projection of 
the wave vector k-- f onto the spacecraft trajectory. 

The same considerations apply to the incident wave at the spacecraft 

- jk • r 

E . = E e (2.62) 

i 0 

CD = k_. * v 

(2.63) 



The signal present at the spacecraft is the narrow-band composite 
of the scattered signals from the entire surface visible from r, with 
each component shifted in frequency by its direction of arrival and 
weighted in amplitude by the root modulus brightness of the surface and 
the antenna gain toward it. Viewed from the spacecraft, the brightest 
region will be the portion of the planet around the specular reflection 
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In addition, the incident wave is also present at the spacecraft with 
the standard Doppler shift. 

It should also be noted that beamwidth, as applied to an antenna, 
and bandwidth, as applied to a filter in the frequency domain, are virtually 
synonymous in this case. The only restriction is that we think of beam- 
width in a plane containing the velocity vector. Since the antenna gain, 

6 ( u )# is 8 function of u alone, and the Doppler frequency for an arbi- 
trary direction is essentially a function of the same vector variable u, 
as k(u • v), it is impossible to have frequencies present from targets 
outside the beam. Conversely, if the input to the spacecraft receiver is 
preceded by a narrow-band filter, it is impossible to have signals present 
from directions such that a) lies outside of the receiver passband. 
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III. INVERSION OF DISCRETE SCATTERERS 


A. HEURISTIC DEVELOPMENT 

In the preceding chapter, we considered the total signal scattered 
from a planet and related this signal to the current distribution on the 
surface. We also found that the magnitude squared of the currents is 
directly related to surface brightness in both an intuitive and rigorous 
sense. Our purpose now is to seek methods of recovering the brightness 
distribution from a measurement of the scattered fields. 

To begin, it is instructive to adopt the standard radar problem 
approach and to consider the radiation from a single point on the planet* s 
surface and the means by which this point might be distinguished from all 
others. This is accomplished by substituting a point source for K(r* ) 
in Eq. 2.38 of the preceding chapter. Since each component of the elec- 
tromagnetic fields propagates in exactly the same way, the vector nota- 
tion for the scattered wave will be dropped. The following equations 
will thus apply to each component of the fields individually and, if 
necessary, the vector fields can always be reconstructed. 

Using two-dimensional Dirac delta function of strength a 

K(r f ) = a 5(r* - p) (3.l) 

-jk. • p - jk — • (r-p) 

1 pr 

E = — (3.2) 


or using the notation of Chapter II 

cp = -k. • P - k — # (r-p) 
r i pr 


(3.3) 


E 

s 


a e 


J<P 


r-p 


(3.4) 
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The geometry is the seme as that given in Fig. 7, with r' replaced by 
p. As a function of time, t, or of the spatial variable s along the 
path, this fits neatly into a line in cylindrical coordinates, such as 
the one illustrated in Fig. 9. 



Here, we have depicted the trajectory of the complex envelope signal 
as the spacecraft approached and moved past our target. Frequency devia- 
tion from the carrier is given by the tightness of the winding with the 
sign of frequency deviation determined by its sense, i.e., right or left. 
Amplitude is indicated by the diameter of the coil. There is an initial 
increase in phase due to the shortening of the path as the spacecraft 
approaches, followed by a stationary point at the moment of closest ap- 
proach, and finally by an unwinding of the coil as the spacecraft recedes 
Suppose that we consider another point p', which is different from 
p. Then, as before, we have 

K(r' ) = a' S(r'-p' ) (3.5) 

and 

a' e JCp ' / \ 

E • = - — - (3.6) 
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This second primed signal is similar in every respect to the first, and 
Fig. 9 represents E' as well as it does E^. However, the difference 
in their phase functions, cp and Cp' , provides an important distinction. 
If, for example, the spacecraft passes the unprimed target before the 
primed, then the signal from the first is rapidly unwinding at the sta- 
tionary point of the second. If the spacecraft passes both targets at 
the same instant so that the stationary points occur together, then either 
the targets are at the same distance from the spacecraft or they are not. 
If they are not at the same distance, then the signals must wind and un- 
wind at different rates since (r - p) ^ (r - p') for any point other 
than the stationary point. Furthermore, these characteristics of the 
phases are independent of the amplitudes 

A method for separating these two signals, i.e., recovering the 
amplitudes of the currents at p and p', suggests itself. Through 
the linearity of the scattering process, the composite signal from p 
and p' is the sum of the contribution from each. 


E 


a e 


jcp 


r-p 



(3.7) 


Therefore, if we can unwind E by the amount corresponding to the 

phase of either E or E ' , the result will be a complex constant plus 
s S 

a term that fluctuates with the difference in phase of the primed and 
unprimed signals. Mathematically, this corresponds to multiplication by 
or e”^' with the result that 

“ j \J/ 

e ~J9 E = _L— + a ' (3.8) 

l r-p | |r-p'| 


e -jCp' E = + (3.9) 

| r-p l | r-p* ( 
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where 


t = <p - <p* (3.10) 

Integration over many fluctuations of e ±3 * should then reduce the 
unwanted terms to a small fraction of the total. For example, from 
Eq. 3.8 we would expect 





FIG. 10. DIAGRAM FOR PROCESSING SIGNAL FROM 
TARGET AT p . 

From this, we may immediately generalize that for a set of discrete 
scatters 
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(3.12) 


E 



a . e T i 
1 



and 


xk 


= CPi - 


9 ,, 


(3.13) 


To recover 



e _ jcpm Edi 




a. 



= am 


(s r s 0 ) ^ 

+ ) (small terms) 

1-pJ ^ 


(3.14) 


These operations are possible because of our assumed knowledge of the 
spacecraft trajectory, 

A block diagram for these operations is given below in Fig. 11. 



FIG. 11. PARALLEL PROCESSOR FOR 
MANY POINTS. 


It should always be borne in mind that E represents data, a mea 
surement of the scattered wave, over some portion of the spacecraft 
trajectory . 
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B. THE RESOLUTION CELL 

The ultimate usefulness of an image lies in the finest detail that 
it preserves, or, in turn, upon the minimum distance that can exist be- 
tween two points and still allow them to be distinguished from one another. 
The classical optical function used for this purpose is the spatial dis- 
tribution in the image space of the field intensity corresponding to an 
optical point source, or the image point response. 

A similar approach is applied to the monostatic radar problem, in 
that it is the radar system’s power response to a point target that is 
taken as a measure of the system’s capability to resolve closely spaced 
targets. Instead of an optical system with an image plane, we now have 
a receiver processing the returning signals in time-frequency space to 
determine target range and velocity. The resolution that may be obtained 
is directly related to the waveform used, and is characterized by the 
ambiguity function. Unfortunately, the ambiguity function does not apply 
here for two reasons: (l) we are interested in a resolution cell on the 

surface of our planet rather than in time-frequency, and ( 2 ) the signal 
arriving from each of our scatterers is a different function rather than 
a standard signal with a continuum of time-frequency displacements, as 
required by the assumptions in the development of the ambiguity function. 

However, the physical basis for calculating the point response is the 
same for both the radar and optical problems. That basis is the assump- 
tion that the absolute phases of all parts of the target are statistically 
independent, so that on the average, it is the powers, or intensities, of 
the responses from each target point that will be superimposed in the 
image. This same assumption will be made here, but a hybrid approach is 
taken to the calculation. The signal is considered as processed in a 
manner similar to that described in Section A above. But, for purposes 
of analysis, it is associated with the geometry that existed when the 
signal was received. With this approach, our problem becomes somewhat 
geometrical . 

The points p and p’ and a fraction of the spacecraft’s trajectory 
are shown in Fig. 12. For the purpose of this computation, the trajectory 
is assumed to be linear. 
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RELATION OF p AND 



FIG. 12. 
P f . 


A coordinate system is centered on the point p as shown. For 
convenience p, p* and the trajectory are also assumed to be co-planar. 
The functions r and r f are the distances to p and p* from the 
point s respectively, while R is the distance of closest approach. 

Any results obtained from the co-planar assumption may be applied 
to the planetary case by rotating points on the sphere about the trajec- 
tory until they lie in a common plane. Using our new variables, the 
analytic signal on s is 


E(s ) 


a e 


J9 



(3.15) 


cp = -kr 


cp' = -kr' 


(3.16) 


Applying the operations of Fig. 10, and taking the magnitude squared to 
obtain power 
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Expanding Eq, 3.17 



2 


a + 



ds 

r 


2 


+ 


2R£ 


<a 


a 1 



(3.18) 


The last term is taken equal to zero by the independent phases argument. 
Thus, we are left with 





(3.19) 


where the integral factor describes the effects of the target at p' on 
our evaluation of the brightness at p. This factor is also precisely 
the image point response, since it is the expression one would obtain 
for the response at p', due to a source at p. Therefore, we are led 
to a calculation of 


: u„ . r f 


(3.20) 


Starting with the exponent 


t = -k(r-r ' ) 


(3.21) 


/2 2 s 2 

r = (R + s ) 


= R 1 + 


2R 


f <<:i 


(3.22) 
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r' = [(R-Tl) 2 + (s-|) 2 ] 


= (R-n) 


i + ^ 


2(R-n) J 


| ^ « i 

R-Tl 


(3.23) 


the difference between these two terms is 


r ' r ' = Is ■ Ifs^T + 1 


_ 2R£s - s 2 T] + R£ 2 + 2Rn(R-n) 
2R(R-n) 


(3.24) 


When applied to the integral, the last two terms in Eq. 3.24 will give 
only a constant phase factor and may be dropped leaving 


r-r 1 = 


-s 2 ri + 2Rs£ 
- 2R(R-:fr 


(3.25) 


as the important part of the path difference. Completing the square in 
Eq. 3.24 yields 


2R(R-r| ) i 


{r\ /z s - rf^gR) 


2 A 2 


T) ^ 0 


(3.26) 


and dropping the constant term again and substituting 2 tt/A for k fur 
ther reduces the exponent expression to 


y(s) 


-n/A 

R(R-n) 


(V 


rf^R) 


(3.27) 


Under the assumptions required for the binomial expansion approxi- 
mation in the exponent, the variation in the denominator of the upper 
integral is also small and approximately equal to R. The lower integral 
may be evaluated in terms of the inverse tangent. 
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~ K / 


•* <*• - **»> 


ds 


R tan 1 s/r 


1 a. 


(3.28) 


letting 


;(s) 


(ar(r-ti)) 


54 u 

(n /2 s - n /e |r) 


(3.29) 


yields 


.<M> — (^f r (Sl) e' J * ’’ a, (3.30) 

“1 /— I 1 V / z 


R tan s/R 


2n ' 


If we now express R tan 1 s/r 
then 


S 1 

as approximately equal' to s - s , 
s 2 10 


g(S» T l) = 


(A Rkll iy y* Z ( s 1 ) 

(vO ' ( COS 2 


z(s 0 ) 


2 ^ , . t 2 \ 

z + j sm ~ z ) dz 


(3.31) 


Thus 


i(l, i) = |gU,n )| 2 = . |[c(«) 

2 t](s 1 -s 0 ) 


= ( Sl ) 


z ( s 0 ) 





z(s o } 


(3.32) 
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where 


and 



. it 2 

sin - z 


dz 


Jt 2 

cos — z 


dz 


are the Fresnel integrals for sine and cosine respectively. 

The expressions may be evaluated on the £ axis by returning to 
Eq. 3.25: substituting 


T) = 0 , r-r ' = (■ s/R 


(3.33) 


SO 


g(S,o) = 


s r s o 


f 




3 /r 


ds 


(3.34) 


Equation 3.34 may be interpreted as the Fourier transform of a single 
pulse starting from s^ and lasting until s^ so 

g(|,0) = sinc^/A • 1 R -J (3.35) 

l(|,0) = sinc 2 ^|/A • 1 R (3.36) 

6 5 

A contour plot of l(|,t|) for R = 10 , = 10 , is given 

in Fig. 13. 
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g(0,Tl) = J J 


f/2 e" ^ f dz 


(3.40) 


Using the Fresnel integral again yields 


<°^ = ^ ( c WH 2 [*(i)]) 


(3.41) 


It is a simple matter to show that this is well behaved at the origin 
and that l(o,0) = 1. The £-axis resolution may also be expressed in 
terms of the quantity f by substitution in Eq.‘ 3.36: 


1(1,0) = sine (f|/A) 


(3.42) 


With these expressions, it is possible to evaluate the resolution 
in the £ and q directions. A natural choice in the £ direction 
is the distance between the first zeros of the sine function. Thus: 


(a|/a) f = 


or 


ai/a = j 


(3.43) 


The choice is not quite so clear in the q direction since l(o,t|) 

2 

has no zeros, and we take the value l(o,q) = 0.1 ~ l/n arbitrarily. 
Then using the asymptotic value of the Fresnel integrals 


C(z) - l/2 
s(z) - l/2 


as z 


gives 


l(0,n) - 2 ~ 2 At] 

jc f 


(I) * (;)' 


A 
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or 



(3.44) 


Doubling the result to account for both the T] and -t\ directions, 


Aa _ *1 
A .2 


(3.45) 


The resolution in the two dimensions £ and T] has been plotted in 
Fig. 14 from Eq . 3.43 and Eq. 3.44. 



FIG. 14. RESOLUTION VS f 
NUMBER. 


Obviously, the key to resolution in the r\ direction is to make 
f as large as possible. Note that so long as our approximations are 
satisfied, the results are independent of range. 

A critical factor in the validity of the point response function 
l(|,T]) as an estimate of the resolution is its behavior for large £ 
and T]. If l(£,T)) does not fall off rapidly enough, then it is possible 
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for a collection of targets far removed from the origin of the 
coordinate system to contribute a significant fraction of the total power 
when the data are processed to recover the signal from a particular p. 
The importance of this effect may be estimated by imagining that the con- 
tributing scatterers are uniformly distributed on the target surface, 
and by calculating the fraction of the total response that comes from a 
small area about the target point. Making the change of variables 


£ = p cos 9 


T] = p sin 0 


(3.46) 


the pertinent quantity becomes 


A(p o ) = 



(3.47) 


It does not seem possible to evaluate the integral from the expression 
for I ( p , 9 ) which we already have, so we revert to the geometrical repre- 
sentation and begin again. 



FIG. 15. GEOMETRY FOR CALCU- 
LATION OF A(p). 
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From Fig. 15, the phase difference is 

t 

V - ^ * p cos(e-ct) 


(3.48) 


and the pertinent integral is 


g(p,e) = ~ 




da, a = 


£ 

R 


ot m a. 


i(p,e) = g • g* = f f “ ^( ar )) ^ da t (3.49) 

*7V *VV 


0 0 


The integral over p, 9 is 


r p o r 2rt r p o /*2n /* a i r a i 

I I i(p. 0) pdpde = I III 

J o o *'o •'a a 


pe 


-J ( v(ce) - 'Ka' )) 


dada'dpde 


(3.50) 


Now 


t(a) - t(a' ) = — p[cos(e-a) - cosCe-a' )] 

A 


(3.51) 


while 


jz cos 0 = J^( z ) + 2 ^ (-) kj 2k ( z ) cos ( 2 k 0) 

k=l 


+ j 2 ^ ( ~ ) kj 2k+l sin [ 2 ( k +l )0 ] (3.52) 

k=l 
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and therefore 


-j p(cos(e-a) - cos(e-a l ) -j 4^ p cos(e-a) + j p cos(e-a r ) 
e A = e • e 


= J o (a 1 p ) + X (sines 


and cosines products) 


(3.53) 


Integrating 0 over 2jx, an approximate bound is given by 


c w ° r 2 r 0 f 1 r 1 2/2* \ 

I p i(p,e) dpae < 2«l J / p j 0 (t p ) dp 

0 J o J o a J a 


da da' 


= 2fl(a 1 o: 1 ) 2 J pJ Q (^ p) dp 


(3.54) 


The standard asymptotic form for J q (z) is 


/ 2 


( 2 \ 

— ) sin z + 0 
jtz / 


z as z -► 00 


J 2 (z) -(-2-) si 
o \Ttz / 


2 

sin z 


(3.55) 


so that Eq. 3,54 is approximated by 


/ r pi(p,e) dpd <-| (a -a 0 ) 2 (p Q ) (3 * 56) 

J 0 J 0 rt 

For an infinite target plane the integrated response does not converge. 

Equation 3.56 may be applied to a finite surface by placing an upper 
boundary on the integral in the denominator of P max ' Supposing that 
the target surface is made up of a set of discrete scatterers of average 
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power P and occupying area A, we have for the power from the total 

area included in o 

K max 


\ <>max (P / A > < 3 - 37 > 

Tt 

n 

The power from a target at p is p(a -a Q ) /2 so that the ratio of the 
two is 


(Power from target) 
(Power from total area) 


Ap 


P (« l^ 0 ) 2 /2 


max 


p/a 


An 


2Ap 


max 


(3.58) 


Note that this result is independent of f number and hence independent 

of the length of the path over which the data are taken and of the shape 

2 

of the resolution cell. The area of integration jto may be limited 

'max 

either by the spacecraft antenna or the size of the planet itself. 

C, THE EFFECT OF ERRORS ON THE INVERSION 

In the last section, we considered the highly idealized case of an 
exact inversion with a linearized trajectory. By exact, we mean that 
our knowledge of the relative positions of the spacecraft and target is 
so complete that we can compute the 9-: without significant error. For 

J 

the calculations made so far, we have, in effect, assumed knowledge of 
the spacecraft trajectory to a minute fraction of a wavelength. At best 
this is a highly unrealistic assumption. Furthermore, if high resolution 
is to be achieved, particularly in the dimension normal to the trajectory, 
it is necessary to have trajectory information over a path which is a 
significant decimal fraction of the range, R, to the target. In prac- 
tical terms, this requires the collection of data over paths that cover 
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hundreds of thousands or even millions of wavelengths. Therefore, we 
consider the effects of errors on the reconstruction to determine the 
uncertainty that can be tolerated in the inversion. 

The motion of a spacecraft is extremely smooth. In the two-body 
problem, where we need consider only radial accelerations, the derivatives 
of the spacecraft motion diminish rapidly with increasing order. Even 
when three bodies are concerned, as, for example, in the effect that a 
local moon or the Sun might have on a planetary orbiter, the effects of 
the third body are perturbations acting over periods of days or weeks. 

(cf. Ref. 17). Consequently, we need only consider the effects of an 
erroneous orbital determination. In other words, the errors in one cal- 
culation of the cp.(s) are due solely to an inaccurate determination in 
J 

our theory of orbital motion. In Fig. 16, suppose that the true trajec- 
tory is s, but we think that it is s'. 



FIG. 16. DEFINING THE ERROR 
FUNCTION FOR THE SPACECRAFT 
TRAJECTORY . 


If we process for the point at p, then we are in error by e(s) rad/m 
and the point response function is 



(3.59) 
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Using the previous approximations for ^(s) 


I e (i,n) a 


s 

/ 


2 

^ rt/A(s T1-2R s) . / \ 
, J R(R-n) — + je(s) 


ds 


(3. 60) 


Supposing that the error is of the form 

e(s) = as 2 + bs + c (3.6l) 

gives, from our expression for 


t + € = 



(3.62) 


Evidently, the result of errors up through the second order is to cause 
a shift in the response function by the amounts Af, Atj 


An = aR(R-n) - 


(3.63) 


That is, nothing is destroyed by errors of this type. The response func- 
tion is preserved, but the signal we associate with p at the origin of 
the |t) coordinate system really belongs with A|, At) . 

Considering only one term in e(s) at a time reveals the physical 
significance of these terms. Denoting the angular error between s and 
s' by u 


b 


= Uk = 

A 


(3.64) 
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For small errors the displaced maxima will still be near the |,T 
origin and rj will be small compared with the range, so that 


A| = u T 1 * ' ( R-T1 ) = UR 

A 


A 


(3.65) 


Linear phase errors therefore shift the image point by angle .O. 

For the second-order term alone, with the same approximation 
R » T], we have 



R 2 a 


( 3 . 66 ) 


so that the quadratic error translates directly into range. In terms of 
the radius of curvature of the orbit, R q , it is a simple matter to 
show that 



where T is the error in our estimate of R q . Thus 


At) = -2 



o 


.A ) 

<>„ + r)2 / 


(3.68) 


In general, one would expect errors from this source to be quite small. 

The constant term makes no difference at all. 

There is one other error of a sort similar to the ones we have just 
considered, and that is an error in the determination of the carrier 
frequency. In this case, the calculated phase paths that we use in our 
processor will all be in error by an amount proportional to the error 
in A. If we let the true wavelength be A Q and the erroneous wave 
length used in the processing be A £ , then, in the calculation of the 
point response, we must use 
SEL-67-042 50 


(3.69) 




r* 




2(r-t)/ 


where we have again employed the approximation of Eqs. 3.22 and 23 
and the geometry of Fig. 12. 

Simplifying, 


♦ 


2« 

A_ 


R + 


2 A / 

2r - r r 11 + 




2(R 




A 


2R 


A £ 2 (r-T) ) 


A, 


(3.70) 


Again dropping the terms that are constant in s, 




2it 

A. 


s 2 (A o /A € )(R-'H) - Rs 2 + 2R£s 


2R (R-t) ) A £ /A o 



r 

o 

/ K\ 1 


2rt . 

S 

1 

iTj & 

\r< 

1 

r-t 

*i 

+ 2R j* s 

A 1 

£ 


2R(R-T] ) 

j 


(3.71) 


The principal effect is to shift change l(|,T]) to I [ l, T l - R(l - A e /A o ) ] . 

The errors that we have considered thus far are peculiar, in that 
they produce systematic effects on our calculation. This is because 
they enter the calculation only once as the parameter of some function 
as, for example, in the use of an erroneous wavelength. A measurement 
error in the determination of the array of the conical elements that 
describe the path of an orbiter would be another. In this second case 
the mistake can be shown to propagate through the orbital computations 
to produce the linear and second-order effects that we have already 
discussed. We have seen that these errors have little or no practical 
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importance. However, the errors associated with the data itself may be 
critical. These may be imagined to arise in several ways, the most obvi- 
ous of which include oscillator instabilities and sampling errors. More 
subtle difficulties would include the effects of any atmosphere through 
which our signals must pass and even the effects of the interplanetary 
medium on the propagation of the illuminating wave. It is supposed that 
all these sources of error will impose an additional random phase modu- 
lation on our data. The phase error function, e(s), thus becomes a 
random variable, which is assumed to be stationary and ergodic . The error 
function may also be assumed to be zero mean without loss of generality 
since a constant may be added or subtracted from the phase without affect- 
ing our final results. Then, as before 


and 


r S i -j(t+e) 

g(S,T]) J ; ds (3.72) 

S 0 


where 


<g(M)> = 




<e“ je >ds 


(3.73) 


< 



) p(e)de 


and p(e) is the probability density function of the error function. 
Similarly, 
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<*(M» ■ <l«| 2 > - / / £ 5 — 4 J ) H > 


ds ds 1 


s o s o 


■ f / 


1 -J(*(s)-*'(s)) R(ssl) ds ds , 


(3.74) 


where r(s-s') = autocorrelation function of e S ^. 

By assuming that p(e) is symmetrical, we may compute these expec- 
tations for g(0,0) and l(0,0) where ^(s) = 0. For tractability we 
take e(s) as being made up of n independent parts. That is, p(e,€') 
is such that r(s-s') = 0 if s-s' is sufficiently large. Then we write 



k=l 


x + jy 


(3.75) 


The factor Ji is taken as the distance required for r(s-s') to 
go to zero. R is still range and a is the strength of the scatterer. 
Defining: 


a = <x> 

s^ = Var(x) = <^x 2 y (3.76) 

2 

s 2 = Var(y) = <y > 


Clearly, 


<p 2 > = <x 2 >+<y 2 > = s x + s 2 + a 2 


(3.77) 
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and 


For large n, 
density of p 


Var(p) = <p 2 >-<p> 2 = s 1 + s 2 


(3.78) 


1 8 

we may use the results of Beckman to write the probability 


p(p) = 


a 

2s, 


S 1 + S 2 2 
4s 1 S 2 P 




X 


00 


I 

m=0 


(-) 


m 


£ I 
m m 



(3.79) 


where I is the modified Bessel function of order m and 
m 

1 , m = 0 

S = « 
m 

2 , m / 0 

The distribution of e determines everything. If € is narrowly 
distributed then Ot - 1 and e has very little effect. If, however, e 
is distributed over a large part of 2jt, then the mean Ot is small, 
while the variances become large. In particular, if the distribution of 
e is uniform over 2rt, or any multiple of 2k, then p(p) becomes 
Rayleigh. For a Gaussian distribution, p(p) becomes Rayleigh as the 
standard deviation approaches a full cycle. When this occurs, everything 
is lost, since the phase information upon which this method depends is 
completely destroyed. 


(3.80) 


D. INTRODUCTION OF THE REFERENCE WAVE 

If one considers the results of the previous section in practical 
terms, it immediately becomes clear that the requisite coherence for 
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avoiding the detrimental effects of random phase errors is not likely 

* ^ 

to be achieved directly with independent oscillators. Spacecraft veloc- 
ities for planetary flybys or orbiters are typically on the order of a 
few kilometers per second, and altitudes of a few hundred kilometers up- 
ward can be expected. Thus, an f ratio on the order of one-tenth re- 
quires at least seconds, or even tens of seconds at a minimum, to achieve. 
If we consider a direct scheme such as that given in Chapter II, Fig. 1, 
for obtaining the complex signal, then the oscillator used to provide the 
sine and cosine functions and the source oscillator providing the illumi- 
nating wave must be extremely stable. Letting T represent the time 
over which the scattered signal is observed, we require that the spectral 
line of an unmodulated signal have a bandwidth of less than l/x, after 
passing through the receiving system. For a carrier frequency of Q this 

requires circuits with a Q so that 

o 

( 3 . 81 ) 

in both the transmitter and receiver. As a numeric example, consider a 
frequency of 1 gc and an observation time of 10 sec. The minimum Q q 
required would be lO 1 ^. 

However, there is a simple expedient that removes this difficulty; 
that is to use the illuminating wave as it is received directly on board 
the spacecraft as a reference. This is possible because (l) the direct 
signal is generally a great deal stronger than the reflected signal and 
(2) as a result of the spacecraft's motion, the direct signal is usually 
separated in frequency from the reflected signal. After extraction from 
the composite direct/scattered signal present at the spacecraft, the 
direct signal is substituted for the carrier e^*, which we used in 
our description of the analytic signal, and is used to obtain a complex 
signal. This may be simply shown by writing the direct and reflected 
signals as separate parts. 

j(flt - k • r ) 

e = a e 1 + e (t) (3.82) 

o s 
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j (ftt-k • r ) ^ 

Obviously, a e is the direct wave. The e (t) is the 

o s 

analytic signal representation of the scattered wave. The effect of this 

operation is to subtract not only the phase term Qt, normally associated 

with the complex signal, but an additional Hk. • r as well. Thus, the 

th 

signal from the i target becomes 


£ 


i 



+ k, 


r) 


(3.83) 


Since all targets are equally affected, the analysis of the resolution 
cell is undisturbed. We merely subtract the -k^ • r from the expanded 
phase and continue as before. Even the arguments for the effects of 
errors hold virtually unchanged . 

The oscillator stability requirements are relaxed considerably by 
use of the direct wave. Suppose that the transmitter phase jitter is 
given by a random term 0^(t) added to its £}t variation. 

<p = fit - k. • r + 0 (t) 

1 r 6 


While the phase of the reflected signal is 

cp = fit - k- - • (r-p ) + 0 € (t-A) (3.84) 

p i 

The term 

(r-p. ) • u- - - I p-r I 
l p . r 1 1 

l 


gives the total differential time delay along the direct and reflected 
paths, neglecting second order retardation, as illustrated in Fig. 17. 
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FIG. 17. TOTAL DELAY PATH. 



After performing the operations of Fig. 10, the analytic signal from a 
point target becomes 


E . 
1 


a i + k t * r + ^ e (t) - ^ e (t-At)) 

i^J e 


(3.85) 


To the extent that 0 (t) and ^(t-At) are correlated, the variations 
that they represent cancel out. Furthermore, this is independent of the 
length of the observation time T. Again, requiring that Q q be large 
enough to guarantee the needed stability, 


Q >-«“ 

o — 2jt c 


(r-p i ) 


u- - 

P i r 


- r-p 


A 


(3.86) 


Therefore, the oscillator stability now goes directly as the path length 
difference in wavelengths. The distance associated with this Q q , C/At 
is referred to as the coherence length of the oscillator. 


E. THE OUTPUT SIGNAL -TO -NOISE RATIO 

Assuming that our radar is operating within the restrictions for 
preserving phase coherence, all the errors considered so far have repre 
sented imperfections in our knowledge of the experimental parameters 
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under which our data are taken, with the conclusion that the principal 

> 

effects can be summed up as a distortion of the coordinate system to 
which the recovered data are referred. However, the most severe problem 
to be expected in any form of planetary radar is that of overcoming the 
receiver noise with signals that have propagated over extreme ranges 
and are received only after reflection from the surface of a planet. 
Therefore, the scattered signal exists in the presence of additive ther- 
mal noise, and the processor must act on signal and noise simultaneously. 
This is the situation we shall now consider for the point target. For 
convenience in describing the noise, the signals will be considered in 
the time domain. 

The method of recovery has been to multiply by the complex conjugate 
of the expected signal from a given target and to integrate. Thus, if, 
as before, the composite signal is E(t), then the processor may be 
written as a convolution with the expected signal from each portion of 
the surface, sampled at the appropriate time 




T=0 


(3.87) 


where OL is an arbitrary gain factor. This may be modeled as a bank of 
linear filters with impulse response E*(-t), -t^ < T < -t^ sampled at 
T = 0. 



FIG. 18. PARALLEL FILTER 
MODEL FOR PROCESSOR. 


Now, to E(t) we add white noise with complex envelope N(t) so 
that for a single filter, our model becomes 
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£(♦) 


*<£> 


I N(t) 

FIG. 19. ADDITION OF NOISE TO PROCESSOR. 

If we let G^(cjd) be the transfer function of the filter with im- 
pulse response E.(t) and assume that the noise is a stationary random 

1 2 2 
process, the output noise power density is N Q / 2 \ G ± \ Q£ , and the out- 
put noise power is just the density integrated over frequency, 

2 N /-+oo 

Noise Power out = a — J | G (co ) | ckn (3.88) 

•^-00 

By the convolution theorem, the filter response to E^(t) itself 
may be expressed as 



a 



(3.89) 


when T = 0, the ratio of the output powers is 


S 

N 



(3.90) 


Recognizing that the integral in the numerator of the right-hand 
side is just the energy contained in the signal, 


59 


SEL-67-042 



E 


dco 


(3.91) 


-+00 

= | l^M 

OO 


gives 


S 2E 
N “ N 

o 


(3.92) 


Filters with an impulse response which is proportional to a re- 
versed version of the complex conjugate function of a given signal (or 

equivalently, have a transfer function which is the complex conjugate), 

19 

are said to be matched to that signal. This definition certainly 
applies to what has been done here. A well-known property of such fil- 
ters is that they maximize the output signal-to-noise in a mean square 
sense, and that the signal-to-noise ratio depends only on the energy 
in the signal and the spectral density of the noise as given by Eq. 3.92. 
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iv . THREE PHYSICAL INTERPRETAT IONS 


A . SYNTHETIC ARRAY 

Suppose that instead of a spacecraft recording signals along its 
trajectory, we consider the construction of an array of infinitesimal 
antennas, co-linear with the trajectory, as in Fig. 20. 



FIG. 20. THE ARRAY ANALOGY. 


The scattered fields arriving at each element of our array would be the 
same as before and given by: 


e(t,n) 


-j [fit-kr(n)] 

a e 

r(n) 


(4.1) 


In Eq. 4.1, we suppose that all the elements are equally spaced with 

distance A ,£, and r(n) is the distance from the target to the n** 1 

element. If 0 is the phase shift associated with the transmission 
n th 

line running from the n element to the summing point and the elements 
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have equal gain, then the output of the array is 


e 

o 


jftt 




-j[k r(n)+0 ] 
n 

1 TO 


(4.2) 


The magnitude of has a maximum if the ^ are chosen so that 

kr(n) + 0 = constant for all n 

' n 


or 

0 = -k r (n ) + c 

r n 

Substituting this result in Eq. 4.2 we obtain 

-j[kr(n)]-kr(n)] 

S e — 

l 



(4.3) 


V i ~ 

= a 2 ^ ‘ 

-j ( n ) 

However, the multiplication by e and the summation of the con- 

tributions along the path are precisely the operations we introduced in 
Chapter III for the recovery of the |a| . The phase delay introduced by 
the connecting transmission lines corresponds to the multiplication by 
e jkr (n ) ^ while the connection at the summing point performs the integra- 
tion. The difference is that whereas an array simultaneously samples 
and sums the incoming fields along the entire path, the spacecraft can 
only sequentially determine the amplitude and phase field at each point 
and in some way record the result so that the weighting by the proper 
phase factor and integration may be done later. In terms of this analogy, 
the use of the reference wave in Chapter II is just a device for calibra- 
ting the phase measurement. 


N 


if R - r(n) 


(4.4) 
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The usefulness of this interpretation is limited by the parameters 
of our analysis. Since we are always interested in path lengths that 
are a significant fraction of the range at the point of closest approach, 
we must always deal with spherical wavefronts from our targets. In fact, 
it is the different curvature of the wavefronts from various ranges that 
gives us resolution in depth. Thus, our targets always lie in the near 
field of the array. Consequently, one of the more powerful tools of an- 
tenna theory, the Fourier transform relations between the fields in the 
aperture of an antenna and its radiation pattern, cannot be applied 
directly to our problem. Still, the idea that the processing represents 
an array focused at close range is a useful concept. 9 

B. THE TRACKING FILTER 

We have previously considered the inversion as being made up of two 
distinct operations, one where the data are collected by heterodyning the 
scattered signals with the incident wave, and a second operation where 
we correct the data for the phase changes from a given target as a func- 
tion of position and integrate the result. For purposes of calculation 
and data interpretation, the geometrical approach adopted in Chapter III 
proves itself most powerful. 

However, perhaps the simplest approach to understanding the process- 
ing for a point target in familiar terms is derived from a consideration 
of the time-frequency characteristics of the signal scattered from a 
point as the spacecraft goes by, and the changes in the processor for 
the same time period. This is the interpretation we shall now explore. 

Re-drawing a simplified block diagram for processing the data from 
a single point with the signals as explicit functions of time, we have 


09 

R 


s mm T 

jkr(t) 


PROCESSOR 


<8> 


r 1 i 

y 

'to 

x. 

1 

# jV«»> 1 

j(kr(t)-k| * 

• 

i(t) ) 

R 


& 


la I 


FIG. 21. SIMPLIFIED PROCESSOR. 
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ik ± ‘s(t) 

Recognizing that the reference signal e functions only as 

a means of providing a coherent reference signal necessary for practical 
reasons, we may give a mathematically equivalent diagram with only one 
frequency multiplication. 


-jkr(t) 



FIG. 23. TRACKING FILTER FORMULATION. 


The overall effect is to extract continuously the Doppler shift associated 
with a particular target and to integrate the result for a time 
In the frequency domain, this corresponds to tracking continuously the 
target signal with a narrow-band filter whose passband is on the order of 

i/tvV Hz - 

It is a simple matter to verify our signal-to-noise ratio calcula- 
tions. The signal power out is |a|^ while the noise power is kTB/2, 
where B = l/t^-t^. Therefore 
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i 


s la 


2E 

N " kT 1 ' N 

2 T ° 


(4.5) 


as before. 

This approach also provides some additional insight into the resolu- 
tion mechanism, particularly resolution in range. Using our straight 
line approximation to the flyby: 


r = {R 2 + x 2 (t)) /z 


(4.6) 


r = 


xv 


ill 


Multiplying by 2jt/A, we have the instantaneous frequency shift of 


m = -2n - 


i/R 


(l + (x/R) ) 


Vz 


(4.7) 


For a single target this gives the S-shaped curve plotted in Fig. 24. 
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As we have already seen, the recovery of the signal that follows this 
curve is accomplished by continuously tracking this signal with a narrow- 
band filter over some fraction of the trajectory. Signals from other 
targets are rejected by this mechanism. Suppose that there are just two 
targets present, neither of which is at the origin of our coordinate sys 
tern. For convenience, take one as displaced parallel to our trajectory 
and the other as displaced in a direction normal to the trajectory. The 
trajectories of these two points and filter tracking the origin are all 
given in Fig. 25. Clearly the resolution parallel to the path can be 



JL 

R 


FIG. 25. TIME -FREQUENCY TRACKS AND THE TRACKING 
FILTER. 

quite good, since only a slight displacement in the £ direction is 
required in order that the signal avoid the filter passband altogether. 
In the T]-direction , however, the signal always passes through the fil- 
ter, so there will always be some contribution to the output from this 
target, although its relative strength may become quite small when the 
filter tracks the desired signal for a longer and longer time. 

One other important property of the method is revealed by this for- 
mulation. As x/R becomes larger, 

(4.8) 

A 

for all targets near the origin. Thus, if the filter is made to track a 
given signal too close to its asymptotic values, then resolution is lost 
as the signals from other targets approach the same asymptote. 
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C . THE HOLOGRAM ANALOGY AND OPTICAL RECONSTRUCTION 

The analogs of the antenna array and of the tracking filter demon- 
strate a striking connection between spatial and temporal signals, i.e., 
an array integrates in space in a manner that is completely equivalent 
to a tracking filter in the time domain. The third analogy is even more 
remarkable, in that it provides an automatic method for performing all 
the operations implied by the first two for all targets simultaneously. 
Furthermore, it is easily explained in terms of the expression with which 
this work really began, Eq. 2.38 of Chapter II: 



r 1 ) e 


-jkl r-r 1 


ds 




r-r 


(4.9) 


We wish to apply this expression to diffraction phenomena that occur 

as the result of a plane wave propagating through a two-dimensional medium 

of varying transmission. We assume that there are no spatial variations 

-jki .r 

in phase associated with the medium. That is, in Fig. 26,, if e 




FIG. 26. SCATTERING BY A 
TRANSMITTING SCREEN. 


is the incoming plane wave, and T(r) is the transmission coefficient 
of the screen and is real, then immediately after passing through the 
screen, the outgoing wave is 
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(4.10) 


-jk T' 

' E^(r' ) = T(r' ) e 

The integral formulation of the scattered wave may be applied by substi- 
tuting the displacement current associated with E^(r' ) 

K(r') = jo>e E^(r' ) (4.1l) 


for current density. The integral must be evaluated on a surface when 
the currents are known, so we take, as the surface of integration, a 
plane lying just above the diffracting screen. Except for an unimportant 
constant, the field at any point beyond the screen is then given by 
Eq. 4.12 


E 

s 



E i (?') 


-jk|r-r' | -jk^r' 

ds 

r- r ' 


(4.12) 


For the proper choice of r, the integral is identical to the 
Eq. 3.20 we developed for recovery of the amplitude from a point source. 

Suppose that T(r f ) is made to correspond to the data, taken by our 
spacecraft. Since any orbit is always a plane curve, this can always be 
accomplished by scaling the orbit onto the diffracting screen, and causing 
the transmission along the scaled orbit to vary with the data. 

Elsewhere, the screen must be opaque. The data must be in a real 
form, such as that given by Eq. 4.13, for a point source originally at 

p. 


T(r f ) = V 


+ j (k | r 1 -p | +k. • r' ) - j ( k | r-p | -k. • r 1 ) 


I r 1 — p I 


| r T -p 


% 


(4.13) 


If the wave number of the illuminating wave is now scaled by the 
same factor and given a direction that corresponds to the negative of 
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FIG. 27. TRANSMITTING SCREEN 
WITH SPACECRAFT DATA. 


the original reference wave, Eq. 4.12 with T(r f ) given by Eq. 4.13 
becomes 


e* . 


+j(k| r 1 -p | -k i * r 1 ) -j(k| r-r ! | +k i -r ') 


ds 


scaled 

trajectory 


r'-p 


— T* * 


r-r 


\S 

J Q 


-j(k|r f -p -k *r 1 ) -j(k r-r f +k i *r f ) 


ds II 


scaled 

trajectory 


r'-p 


r-r 


(4.14) 

If r = p, part I of this expression is identical to the starting integral 
for the calculation of the image point response. Part II represents a 
spurious response. The intensity of the diffracted electric field due to 
I exactly recreates r, the image point response we calculated in Chap- 
ter II . 

Furthermore, the response occurs at the point in space that corre- 
sponds to the original target position. We must be a little careful. 

The important portion of the integral is the function r - r* in the 
exponent of the kernel. As was pointed out in the original discussion 
of the image point response, all targets lying on the same circle about 
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the original trajectory produce a response at the same position, since 
they are all related to the trajectory through the same phase function. 

It is only by knowing the location of the target surface that we are able 
to map the surface brightness distribution. 

Now we have the opposite situation. For a straight line trajectory, 
the diffraction screen reproduces a bright spot not only at the original 
point, but at every point on a circle centered on T(r) that contains 
that point. However, if we sample the field intensity on a surface that 
has the same geometrical relationship to the T(r) as our planetary sur- 
face had to the original spacecraft trajectory, we will recover the sur- 
face brightness distribution spatially undistorted, albeit, it is convolved 
with our point response function. It is a simple matter to show that for 
a collection of point scatters, the surface brightness distribution is 
convolved with the point response in the reconstruction. Because our 
functional relationship has remained unchanged, the noise behaves in ex- 
actly the same manner as before and processing gain is retained. 

An interesting point is that when r is far removed from r f , so 
that the vector directions (r-r ! ) and (p-r ! ) are approximately the 
same, the integral becomes 


-jk r-p 


r - P 


L 


(l) ds 


scaled 

trajectory 


(4.15) 


which looks like a plane wave originating at p. This approximation 
always applies to II, since -p represents a point behind, i.e., on the 
side of the illuminating wave, the diffracting screen. Therefore, II 
represents a uniform background intensity that apparently comes from a 
point source behind the slit. An obvious method for avoiding this back- 
ground contribution is to choose the geometry so that the real response, 
due to I, falls outside the wave due to II. 

There is nothing in our derivation thus far that is critical in 
terms of frequency. Herein lies a possibility for analog reconstruction 
of the surface brightness. By scaling the data from radio to optical 
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frequencies, it becomes possible to create a diffracting slit, such as 
we have described, by scanning photograph plates with a beam of light, 
intensity-modulated in this manner automatically generates the diffrac- 
tion pattern we have been discussing. A second photographic plate, placed 
in the original target plane, as in Fig. 28, scaled by the same factor, 
will then automatically record the reconstructed brightness distribution. 



FIG. 28. RECONSTRUCTION VIA 
TRANSMITTING SCREEN. 


A two-dimensional diffracting screen with T(r) corresponding to 

the intensity of the interference pattern generated by a coherent source 

of illumination and a scattered wave is called a hologram. It is easy 

to show that the interference pattern contains, in addition to several 

spurious products, the same data that our spacecraft collects in one 

24 

dimension by product detection. Much of the theory of holograms is 

germane to the reconstruction of the surface brightness distribution and 

u J „ 24,25,26 

has been discussed elsewhere. 
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V. RELATIONSHIP BETWEEN THE RESOLUTION CELL AND SIGNAL -TO -NOISE RATIO 


A. PROCESSOR FOR MAXIMUM-LIKELIHOOD ESTIMATE 

The conventional radar problem rather naturally divides itself into 
two parts: (l) the detection problem, where one is concerned with de- 

ciding when a target signal is present in a background of noise, and 
( 2 ) the resolution problem, where, usually with the assumed absence of 
noise, one is concerned with fundamental limits on the radar’s ability 
to discriminate between closely spaced targets. When stated in this way, 
the radar problem is extremely well understood and we have drawn upon 
ideas of detection and resolution in calculating the s ignal-to-noise 
ratio and image point response. 

A little reflection will show, however, that the ordinary develop- 
ment of the detection and resolution problem cannot be carried much fur- 
ther here. Detection, in the sense of determining the presence or absence 
of a target is completely out of place, because one of the major assump- 
tions of this report is that the target is known to exist and moreover 
that our knowledge of its position involves only a small uncertainty. A 
straightforward application of the resolution cell results is inappro- 
priate because the assumptions of high s ignal-to-noise ratios made in 
most developments and of a knowledge of the target environment are gen- 
erally never true for the applications envisioned. Nevertheless, we shall 
still find the developments of Chapter III extremely valuable later on. 

And so, we take a more abstract approach. This time, we assume that 
the brightness distribution over the surface is a random process generated 
in some unknown manner by the interaction of the incident wave and the 
surface itself. Further, assume that we measure the scattered fields at 
our spacecraft as before, and that the relationship between the scattered 
fields and the current distribution is known and is given by our integral 
expression Eq. 2.38. What, then, is the best estimate of the surface 
brightness distribution averaged over some area, and how does the quality 
and form of the estimate vary with the size of the area itself, the sta- 
tistical parameters of the brightness distribution and the noise? In other 
words, how does one trade resolution for signal-to-noise ratio? 
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Probably the most difficult portion of this problem is the selection 
of a stochastic model for the surface with some basis in reality, but 
which is tractable at the same time. Almost any presumption about the 
surface is likely to be wrong. Therefore, one is led to select some 
model for the surface, derive the expression for an estimate on the basis 
of the model, and then examine the sensitivity of the formulation to the 
original assumptions. This is the goal of this chapter. 

We begin by assuming that the current distribution has the same 
basic character everywhere, but that for one reason or another, varies 
in amplitude from place to place. Such an assumption might be based on 
the belief that the same sorts of material stresses and strains act to 
roughen the entire surface but that the brightness varies from place to 
place because of shadowing or variations in the electromagnetic properties 
of the material. 

Designating the complex surface current distribution as a(g,T)), we 
define the average brightness over the area A by 

A 2 = 7 f (a) 2 dg dr) ( 5 . 1 ) 

A 


and then take 


= a ( £ » "H ) /A 


(5.2) 


so that 


(5 - 3> 

The quantity A then is fixed for any given Introducing the 

geometry given in Fig. 29, the signal received along the trajectory £ 
may be obtained from Eq. 2.38. Designating the kernel by f(Ci Sp 7 !) and 
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FIG. 29. GENERAL GEOMETRY FOR 
EXPANDED RESOLUTION CELL. 


neglecting the constants fCC;!* 11 }) 



and 



!,n) d| dn = A 



r(e,n) t(C; l.n) d| di) 


(5.4) 


so that the received signal is 

S(C) = V(C) + N(0 (5-5) 

where n(£) is the complex envelope of the receiver noise when the space- 
craft is at position The quantity N(£) is taken as stationary, 

white, and Gaussian. 

Autocorrelation will be designated by ft with a prefix to indicate 
the corresponding variable. Thus 

0(6, 6 '. W ) = < 7 ( 6 , *0 7 *(l;V)> = y q>(p) 

so that 

Cp = A 2 cp where p 2 = (£-| f ) 2 + (tI-ti 1 ) 2 (5.6) 

a Y y^ 
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and 


n^e.C') = N to B(C-C) 

V <P(£.C) = <v(C) v*(C')> 
s <p(C.C) = <s(C) s*(£') > = n 9 + v cp 

The autocorrelation function of v(£) is calculated from 
v 9 =<v(Ov(C')> = A 2 /* f y(l, T l) 7 *(|, T])f (C 5 5 ,T) ) f*(C;S,n') d^dt'dr)' 


= A 


JaJa 7^ 


ff d|dt]d| 'dtj 1 


(5.7) 


The term is clearly Hermitian and it is a simple matter to show 

that q> is also. Let g be an eigenfunction of cp so that 

V II. 


/• 


g g = 5 
p v pv 


(5.8) 


and 


A g 

p v 


= j v ^ v ^ 


(5.9) 


The function <p is also positive-definite since 


/h(t) v(t) dt I rff h(t) h*(t')<(v(t) v*(t'))>dt dt' 


= JJh(t) h( t 1 ) v cp(t , t 1 ) dt dt 1 > 0 (5.10) 
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for all h(t ) . 

Therefore, by positive definiteness and Hermitian properties of 

27 28 

^Cp, the are positive and real. It follows from Mercer* s 9 

theorem that the autocorrelation function may be expanded in terms of 
its own eigenfunction as 


v 


cp = 


y a g 


i 


(5.11) 


For 


S 




+ Y 


(5.12) 


we have 


/S S*> = f f g*g <SS*>dt dt' = 2 (a 2 A + N )6 (5.13) 

\ H v / J J (Ji v x x p. o nv 


thus 


<S S > s 0 
N |i v / 


(5.14) 


A linear operation on a Gaussian variable produces another Gaussian 
variable. Thus the S are independent, Gaussian random variables of 
zero mean and 


<K i 2 > - < x !> - < y !> - ‘('V'.i 


(5.15) 


which implies 


pSy* 2 ) =n— r 1 :• 

2 rt (A A + N ) 

H o 


. y l s J 2 

1 + N o> 


(5.16) 
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To obtain a maximum likelihood estimate, we require 




(5.17) 


In turn, 


I 


A [|S | 2 /2 - (A 2 A +N ) ] 

" ° 


(5.18) 


29 

This relation has also been given by Hofstetter for a similar 
estimation problem, starting from our Eq. 5.5. 

2 

In general, Eq. 5.18 cannot be solved for A . However, we make 

2 2 2 

some progress by expressing our estimate of A , A in terms of A 
itself 


or 


I d^l 2 / 2 - 


A = 


(N + A A ) 

_2 k ±_ 


I 


i" < A \ * N o )2 


(5.19) 


I <I S J 2 / 2 - "o’ 


A 


A = 


(N o /A 2 + a m ) 2 


A 


! <V A + V‘ 
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where 


N 

o 

2 

A 


is 


(noise power per unit bandwidth) 
(signal power per unit surface area) 


This will become a critical parameter. 

i i 2 

Substituting for S | from Eq. 


5.13 


A 



00 



1 


(a 2 A +n) 2 

U O 


I 


«*\ * -o) 2 



(5.20) 


so the estimate is unbiased, at least as long as our choice of A in 
the processing is correct. 

While the general variability of the processor form with changing 
2 

values of A is of great interest, the case of special interest to us 
is 


N 


(5.21) 


so that the estimator becomes 


/\ 


A 


2 


I <|Sj 2 /2 - N o ) A^ 



1 


(5.22) 
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For this case, the variance of the estimate is easily shown to be 


* 2 (n o + *V\ 2 . 2 ( v a2 + v ) 2j \ 

Var (A 2 ) = — S ^ = A 


2 1 


[? -:I 


00 t 

[!>;] 


(5.23) 


Hofstetter has computed the efficiency for this case with the result that 


eff = 


2 


A 


N 


2 


a 2 a 


1, -% - » (5.24) 


1 + 


N 


v 2 2 


1 (l +S-* 

V N o >V 


so that for N q /a large, this estimate approaches the Cramer-Rao lower 
bound for the variance of the estimate. 

The estimator is basically 


A 2 = k 


2 


S | 2 A 

pi p 


(5.25) 


Again, for the case of interest, the variance is 


Var A 


N y A 2 

o p 


(5.26) 
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The autocorrelation function is the important ingredient. Maximum- 
likelihood processing consists, at low signal-to-noise ratio, of matching 
the processing filter not to the expected signal, but to the expected 
autocorrelation function. 
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B. CALCULATION OF SOME SIGNAL AUTOCORRELATION FUNCTIONS 

Having determined the general form of the estimator in a fairly 
unrestricted manner, we are now obliged to make some assumptions about 
the surface in order to proceed further. Since, as has already been 
pointed out, the usefulness of any of the results depends on the sensi- 
tivity (or rather insensitivity) of the estimator to our surface model, 
several different cases will be considered. 

The estimator is a time-varying filter whose impulse response is 
equal to the expected autocorrelation function of the scattered signal 
As before 


v <p(C»C') 



*(£;*), T) d£ dr) df 


dTl' 


(5.31) 


where the various terms are defined by Fig. 30 below, and we have assumed 
a linear trajectory once more. 





e 


-jkr 


r 


and 




jkr' 

e 

r' 


(5.32) 
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Considering the most extreme case first, 


we take 


and 


cp = 5(p) , p 2 = (AC 2 + ATI 2 ) 

7 

(5.33) 

AC = I - S' 

At) s t) - t) ' 


This would correspond to a "white" surface. We have 


V <P<C,C)= /Up) 

in ~ rS 


-jk(r-r') 


rr ' 


d| dT) d£' dT)' 


-L 


-jk(r(C;i,T)) - r'(C ; i ) ) 


rr 


d|, dT) 


,<p(C.C’) 


-jk(r(C;|,Tl) - r’(c ' » ^ )) 


As 


rr 1 


as 


^ As 


( 5 . 34 ) 


Thus, as the area goes to zero, we obtain the matched filter used in the 
image point response. To see this 


- X ■■**■/ 


s(C) ,*(C) --(C)) 

r , r ' 


(5.35) 




■(0 ^ kr( . C 2 a ; 


as before. 
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I 


For this case, the variance is trivial 


^.C )! 2 - 


(rr'V 


(5.36) 


so that 


Var A 

(A 2 )' 


/ N \ 

~ I — 1 

\a 2 / 


2 .2 
r r 1 


'N 


//<■> 


dC dC 


JO 

2 , 


«.V 




(5.37) 


Changing changes the estimator. Consider ^ as a linear region of 
width Atj between g = a, £ = b, rj = o. 


„¥«;.«:•) -f r'(C)) dl 

*'6= a 


(5.38) 


Over a short path, the approximations indicated in Fig. 31 are suitable, 


£ 



FIG. 31. APPROXIMATION FOR 
EVALUATING CORRELATION 
FUNCTIONS . 


with the result that the path difference is given by 


£ 



<<-«*> I 


(5.39) 


Therefore, 

i 
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/-*1 ijk((C“?' )^/ R + r (C)-r'(C' )) 
v q)(C,C) = AtlJ * <*fi 

^ n 


ATJ e -Jk (r( C ) - r ' ( C ,)) f S 1 e +Jk(C-C')S/R d6 

h 


-jk(r(C)-r'(C')) 


= Ari 


sine v 


rr 


(C“C l ) 


> A 


where 


2 = £ — E 
5 i 


(5.40) 


A similar extension applies in the tj direction. Having integrated 
over one strip, we need only add up the strips. Again obtaining the 
approximations from Fig. 31. 


v »(c.c> 


.-Jt(r-r') >Jk ( ? . ?; ) e/B e Jk(v/?V 



d£ dr] 


where 


-jk(r-r') /k£ 

= e sine l — 



= t^-tj 


0 



(5.41) 


The appearance of the sine function as a weight on the point processor 
is interpreted as accounting for beats between different portions of the 
resolution cell when it is expanded. 

But how does the processor change with the autocorrelation function? 


y 


9 


-alS-S'l ~ a | r| -t) ' 
e 


(5.42) 


Then 





.-“M - *M dn ^ Bfl, 


(5.43) 
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as before 




- r* 1 


T - C 


(5.44) 


giving 


Cp = H H 


-jk(r-r f ) 4 

e — sme 

a 

_*f - VP^Jsincff («■)] 

.* ♦ ft ■ )\ 




provided a is large compared with the area of integration, 
the simplest case to visualize is the Gaussian 


(5.45) 
Probably 


9 = e 


- f(A £ 2 + An 2 ) 


which yields 


< 2 + c * 2 


cp = HH 


e 


a 


kH „2 „2 

2 Slnc C +R 


C' 2 +R 2 sin[^(C-C')] 


(5.46) 


Finally, the coherent surface gives 

9=1 

1 


Cp =H 2 H 2 Sine ^ (Je^ - ^ ,2 +R 2 ) sinc [jf (C“C')j (5.47) 


C. THE SIGNAL -TO -NOISE RATIO 

Under the assumption that thermal noise dominates the process made 

A A 

9 2 

in the derivation of the variance of A 45 , the ratio of the power in A 
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is the signal-to-noise ratio. That is, 


/\ 

to the variance of A^ 



2 N J J 
Var A o 


(5.48) 


The autocorrelation function cp calculated in section B can be used 

v T 

in the integral. For convenience, the expression will be normalized by 
2 , 4 

N /A and only the integral will be determined. 
o' 

First, for the completely incoherent surface from Eq. 5.33 


cp = 6( p ) 
7 





(5.49) 


Substituting for 


//A 4 • S/», 


ffw \ 2 , s v// sl '„c 2 [!f <«.)] al 0 c 2 [f(iw 


dCdC 

(5.50) 


This is easily approximated for expansion along £ by 

JJ |<pK d^ ~ h 2 h 2 (q-y 2 for h, h « (q-c o ) 

~ H 2 H (q -^) 2 for H « (q-^ Q ) , 

H < (C,-^) . 

and H < E (5.5l) 


Similarly, the coherent surface gives 
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<p * 1 

7 


| v cp| 2 = H 4 H 4 sine 2 [f («• )] si„c 2 [f(^? - 


(5.52) 


so that 


//iv'Pl 2 d C ~ H 4 H 4 (C r C 0 ) 2 for H, s « (C r C 0 ) 


h 4 S 3 (C 1 -? 0 ) for H « (Cj-Cq) ■ 

h < (q-y . 

and H < E 


(5.53) 


The Gaussian case is also instructive: 


7 


■p = * A ^> 


I I ^ tt 2 2 e 
I cp| = H * - 


_ £j£L 2 

a/2 


sine 2 (C"C )] sine 2 |p^4; 2 +R 2 " 7c ,2 +R 2 ^ 


(5.54) 


2 

If l/a « H,H, then the sine dominates the integral so that 

//m 2 <1? dc ~ (c-y ( 5 . 55 ) 

This corresponds to a correlation area that is less than the size of the 
area & . When compared with the completely coherent expression, it 
immediately follows that 

2 __ 

Hr (?!-?„) <H 4 H 3 (C 1 -C q ) (5.56) 

a 

as one would expect. 
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VI. REALIZATIONS OF THE ESTIMATOR 


In this chapter we will show how the kernels derived in Chapter V 
may be realized. The form of the maximum-likelihood estimator 


/ 


* * 
scp s 


dC dC 


derived in the last chapter, is equivalent to the optimum detector for 
stochastic signal . ^ ^ In the succeeding few paragraphs we shall 
lean heavily on the detection problem work but shall not follow it ex- 
actly, and we shall put it in a form better suited to our purposes. 

It is convenient to have the form of the estimator in real variables. 
Breaking up the correlation function and signal into real and imaginary 
parts , 


S = t + JTJ 


cp = y + j© 


( 6 . 1 ) 


we can easily 


show that our estimate becomes 


/\ 




[l(t)!(s) f(t,s) + T](t)Ti(s) ¥(t,s) + 2£Tie(t,s)]ds dt 

( 6 . 2 ) 


In obtaining this result, we have made use of the facts that 

1mA = 0 

and © = © since ^ is Hermit ian. 


There are several alternatives available for obtaining A other than 
direct calculation of the integral. Since in the development of cp we 
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“ jk C “i* ^ ) 

always had an e factor, will always contain cos k(r-r' ) 

and 0 the sin k(r-r'). Denoting the other factors in ¥ and 0 
collectively as h(s,t), the kernel is 


l(t) |(s) cos k(r-r') h^(t,s) + rj( t ) t](s) cos k(r-r') h^(t,s) 

+ 2r l(t) £(s) sin k(r-r' ) h (t,s) (6.3) 

u 


Assuming that h^(t,s) = h(t-s), as it generally is from the previous 
section, the realization is 



FIG. 32. FIRST REALIZATION, TIME INVARIANT 
FILTER. 


A second, more general, realization is given by Price and Green 
33 

from Middleton. For this case, we must have an impulse response sat- 
isfying 


h(t-o-) = J 


h 2 (a f , t) h 2 ( cr 1 , ex )d a' for all t in (a,b) ( 6 .4) 


Then we have 
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ijtt) 



FIG. 33. SECOND REALIZATION, TIME -VARYING FILTER. 

The third realization is the familiar weighted radiometer. We require 
a filter h 3 (ir) so that 


/*+oo 

|H 3 (cd) | = I h(cr) e Ja)cr dcr 

■'—oo 


(6 


Then we have 



FIG. 34. THIRD REALIZATION, WEIGHTED RADIOMETER. 
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The fourth realization is in the form of three correlators followed 
by a special filter’. The deviation follows directly from a change of 
variables in the double integral. Taking the first term as an example, 

£(t) £(t') cos k[r(t) - r' ( t * )]h(t-t' ) (6.6) 


Changing variables 


t - t' = A 
t + t ' = & 


(6.7) 


we have 


A = 


t -t 
1 0 


A = - 0 


h(A) di5 dA 


( 6 . 8 ) 


which is achieved by the processor given below. 




FIG. 35. FOURTH REALIZATION, AUTO -CORRELATOR . 
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All of the above realizations may be achieved physically by means of 
regular filtering techniques or numerically with digital computation. 

This latter possibility offers considerable flexibility since the require- 
ment of physical realization may be removed. 

There is still at least one case of some interest left; this is the 
holographic or diffraction screen realization of the maximum-likelihood 
estimator. As an example, consider the situation in Fig. 36. This is 




FIG. 36. FIFTH REALIZATION, OPTICAL ANALOG, 


essentially Fig. 28 for the hologram reconstruction, except that the 
illumination has been modified with a slit so that an extra weighting 
function appears as a factor in the reconstruction integral. Supposing 
that the slit is adjusted to produce the familiar sine (x) diffraction 
pattern, we have for the reconstruction 



sine 


2 


(x-A) dx 


(6.9) 


or 



-jk(r-r' ) 


g(x) g(x' ) sine (x 


-A) sine (x-A) dx dx' 


SEL-67-042 


92 



which will be recognized as equivalent to the Middleton processor for a 
physically unrealizable filter. 33 The integration with respect to A 
may be achieved by moving the slit parallel to the data string at constant 
velocity and recording on photographic film. The film automatically 
squares and integrates. 
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VII. EXPERIMENTAL RESULTS WITH LUNAR ORBITER SPACECRAFT 


A. NATURE OF THE EXPERIMENT 

The first experimental tests of the application of bistatic-radar 
to celestial bodies were made in October 1966, and March 1967, using 
Lunar Orbiters I and III, respectively. These spacecraft carry a te- 
lemetry transmitter which, with appropriate filtering equipment on the 
ground, can be used as a source of continuous wave radiation on a wave- 
length of 13 cm. In addition, these spacecraft each have two antennas, 
a discone omni-directional antenna and a high-gain paraboloid, which 
can be arbitrarily oriented in space. Observations were made using both 
of these. 

The radar parameters for the two cases are given below: 



Omni-Antenna 

High-Gain 

Paraboloid 

Radiated Power 

100-mw CW 

3-w CW 

Antenna gain , transmitting 

-1 db 

22 db 

Antenna gain, receiving 

57 

57 db 

System temperature 

W 

0 

O 

00 

160 °K 


The polarizations of the two transmitting antennas were not deter- 
mined; the polarization of the receiving antenna is right-hand circular 
in the radio sense . 

The receiving system temperature was largely dependent on the phase 
of the moon during the observations and this, along with certain con- 
straints on the Sun, Earth, and spacecraft geometry for the use of the 
high-gain antenna, accounts for the differences in system temperature 
for the two cases . 

B. EXPLANATION OF THE OBSERVATIONS 

Four examples of the records obtained are given in Figs. 37-40. 
These records represent spectra, measured with respect to the instanta- 
neous spacecraft carrier frequency. Time advances to the right. Signal 
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FREQUENCY DEVIATION 
FROM CARRIER (KHz) 


STANFORD UNIVERSITY 

1 



2295 MHz 
CARRIER 


2308 2309 2310 

TIME (GMT) 


OCCULTATION 

IMMERSION 


LUNAR ORBITER OCCULTATION 12 OCT 1966 
100 MW RADIATED POWER 
OMNI ANTENNA 

FIG. 37. TIME FREQUENCY TRACKS, LUNAR ORBITER I, 12 OCTOBER 1966, 23 GMT. 
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FREQUENCY DEVIATION 
FROM CARRIER (kHz) 



APPROXIMATE TIME (GMT) 


LUNAR ORBITER III 19 MARCH 1967 
3W RADIATED POWER 
HIGH GAIN ANTENNA 

FIG. 38. TIME FREQUENCY TRACKS, LUNAR ORBITER III, 19 MARCH 1967, 0300 GMT 
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FREQUENCY DEVIATION 
FROM CARRIER (kHz) 



APPROXIMATE TIME (GMT) 


LUNAR ORBITER m 19 MARCH 1967 
3W RADIATED POWER 
HIGH GAIN ANTENNA 

FIG. 39. TIME FREQUENCY TRACKS, LUNAR ORBITER III, 19 MARCH 1967, 0600 GMT. 
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FREQUENCY DEVIATION 
FROM CARRIER (kHz) 



-15 


0821 0822 0823 

APPROXIMATE TIME (GMT) 

LUNAR ORBITER m 24 MARCH 1967 
3W RADIATED POWER 
HIGH GAIN ANTENNA 

FIG. 40. TIME FREQUENCY TRACKS, LUNAR ORBITER III, 24 MARCH 1967, 0800 GMT. 
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intensity i, represented by relative shades of gray, with the strongest 
signal, being white. The ordinate is frequency difference fro. the 
carrier in kHz. System noise appears as a speckled background, while 
the weak, continuous horizontal traces are due to imperfect filtering. 
The carrier frequency is the strong, continuous, horizontal trace appear- 
mg either near the top or bottom of the records. 

The observations in Fig. 37 were made on 12 October 1966 while Lunar 
was approaching occultation immersion over the lunar crater 
Kastner on the eastern limb. Radiation was through the omni-antenna, 
and the radar parameters given for this condition apply. The altitude 
Of the spacecraft was initially about 60 km, and decreased continuously 
to about 40 km during the period spanned by the record; its velocity was 
approximately 1.9 km/sec away from the Earth. Reflections from the sur- 
face appear as time-varying signals of decreasing frequency, (in this 
and the succeeding records, it is presumed, because of their continuity, 
that the persistent reflections are due to a fixed scattering center.) 

spacecraft carrier disappears when the vehicle moves behind the limb 
at the time indicated as occultation immersion. 

Figures 38 and 39 represent observations of Lunar Orbiter III made 
on two separate orbits on 19 March 1967. Radiation was through the 
high-gain antenna with the appropriate radar parameters. The antenna 
direction was to the side of the spacecraft's flight, and such that it - 
intersected the moon with an angle of incidence of approximately 45°, 
and at a range of 600 km. The geometry was such that the received energy 
corresponds to oblique scattering. As one would expect from the develop- 
ment of Chapter II, the reflected energy is narrowly confined in frequency 
by the limits of the antenna beamwidth. Again, in Fig. 39, several dis- 
crete scattering centers stand out, and their progress through the antenna 
beam is easily followed in the frequency domain. 

The conditions under which the records in Figs. 38 and 39 were made 
differ, in that the beam was directed at slightly different portions of 
the lunar surface. In the first, i.e., Fig. 38, the beam was directed 
at the crater Langrenus at the time the high-gain antenna was switched 
on, and the crater was then slowly carried out of the beam by the motion 
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of the spacecraft . Figure 38 Is also of Interest because the specularly 
reflected signals from the side lobes of the hlgh-galn antenna are visible 
just below the carrier at the start of the record. For Fig. 39 the beam 
was directed to pass just south of Langrenus across an older region of 
the surface, which is composed primarily of maria. 

Figure 40 gives observations under circumstances similar to those 
in Figs. 38 and 39, but the antenna was made to sweep an area of the Moon 
from the crater Albategnius near the center of the disk, south of Delambre, 
near Maskelyne and into the Mare Tranquillitatis . The features pointed 
out in Figs. 38 and 39 are present in this record as well. 


C. REDUCTION OF OBSERVATIONS 


The record in Fig. 37 has been used to demonstrate the nature of the 
two-dimensional information inherent in these observations. The frequency 
differences shown on the ordinate of these records are the instantaneous 
rate of change of phase with respect to the carrier, or direct signal. 

The quantity measured is 


while 


Af - Af R Af D 


Af = — sin 9 
R A 


(7.1) 


The Doppler shift of the direct signal, may be calculated from the 

known orbit of the spacecraft, thus allowing the angle between 

ity vector and the vector to the target, 9, to be determined . Any given 

frequency offset therefore determines a cone about the velocity vector of 
half -angle 9, in which the scatterer must lie. The intersec 
this cone with the surface of the Moon in turn defines a locus of points 
on the surface which might contain the scatter. The intersection of sev- 
eral loci, each determined from separate measurements on a single 
at different times, provides an estimate of the locations of the scatter 
on the lunar surface. It is a simple matter to show that this is equiv 
alent to the use of a linear phase term for determining the azimuthal 
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position of the scatterer and quadratic term in the rate of change of 
phase to locate the scatterer in range, as suggested in Chapter III. 

While this represents an application of the idea presented in Chapter III, 
it is not a complete use of the method, since in this case we are using 
a number of estimates of instantaneous frequency rather than a filter 
matched to the entire record to obtain the results. Thus it has been 
demonstrated that two-dimensional information is present in the scattered 
spectrum, but the full potential of the method has not yet been realized. 

Figure 41 is a drawing of Fig. 37 in which the traces have been 
identified by roman numerals. Figure 42 is a plot of the loci associated 
with trace III in Fig. 42 projected onto the T)| plane of the standard 
orthonormal selenographic coordinate system. 34 For purposes of illustra- 
tion we show the intersection on only one side of the orbiter's ground 
track on the lunar surface. One division in Fig. 42 represents 17 km on 
the Moon. This particular record shows not only the locations of the 
scattering center responsible for the trace, but that the center itself 
is resolved into two different scatterers, one that contributes to the 
initial portion of the trace, and a second that comes into view near the 
end. 

The locations of the scatters are given in selenographic latitude 

and longitude in Fig. 43 by the centers of the squares. The size of the 

squares depicts an estimate of the standard deviation of the intersections 

of the loci. These results have been published elsehwere, along with a 

preliminary analysis of the probable scattering mechanisms responsible 

35 

for the observations. 
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VIII CONCLUSIONS 


In this study we have shown that astronomical bistatic-radars opera- 
ting between a ground terminal and a vehicle in space may be used to probe 
and study planetary surfaces. 

The illumination of a planet with radio waves produces a set of 
scattered fields above the surface, which, when compared directly with 
the incident wave, may be processed to produce a two-dimensional map of 
the radar brightness distribution of the surface of the planet. A crucial 
factor is the use of the indicent waves as a reference when obtaining the 
data on the scattered fields. 

The resolution cell is variable in size and may be expanded to in- 
crease the signal-to-noise ratio at the output of the processor. But, 
for any given resolution cell, the inversions derived represent the 
highest signal-to-noise ratio — in the sense of minimum variance — that 
may be achieved with a linear system. 

For a point source, the processor has analog interpretations as an 
array of antennas, a tracking filter, and a one-dimensional hologram. 

The processor for an expanded resolution cell is realizable in a 
variety of ways, some of which represent physical filters, some of which 
do not. 

The hologram analog may be modified to realize the processor in the 
case of certain special choices of geometry for the expanded resolution 
cell. 

In practical terms, it would seem as though the methods for collect- 
ing and processing radar signals that have been given here can in fact 
achieve most of the advantages claimed for bistatic-radar in the intro- 
duction to this paper. If we consider the uplink case, a spacecraft 
instrumented with a relatively simple, narrow-band, phaselocked receiver, 
an antenna, some sort of data storage device, and telemetry readout 
is sufficient for this experiment. Large transmitters and antenna sys- 
tems to provide the illumination are already in existence on the ground. 

Depending on the choice of radar parameters, such an experiment 
would be expected to yield, in addition to its imaging features, data 
on surface and subsurface structure on the scale of a wavelength. With 
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further development of scattering theory and experiments along lines of 
research already begun, it may even become possible to deduce the com- 
position of a surface from a remote measurement of its scattering proper- 
ties . 

In view of the difficulties associated with the optical mapping of 
Venus, the potential application of the imaging technique to the study 
of the surface of that planet is of special interest. 


SEL-67-042 


106 


REFERENCES 


1. Skolnik, M. I., Introduction to Radar Systems. McGraw-Hill, N.Y., 

1962. 

2. Skolnik, M. I., "An Analysis of Bistatic Radar," IRE Trans., ANE-8 
19-27, March 1961. 

3. Eshleman, V. R., "Radar Astronomy of Solar System Plasmas," SR No. 4, 
Stanford University, 1964. 

4. Fjeldbo, G., "Bistatic-Radar Methods for Studying Planetary Iono- 
spheres and Surfaces," SR No. 2, Stanford University, 1964. 

5. Fjeldbo, G. and V. R. Eshleman, "The Bistatic-Radar-Occultation 
Method for the Study of Planetary Atmosphere," J.G.R. 70, 3217-3225, 
1965. 

6. Kliore, A. J. et al., "Occultation Experiment: Results of the First 

Direct Measurements of Mar's Atmosphere and Ionosphere," Science, 149, 
1243-1248, 1965a. 

7. Fjeldbo, G. , W. Fjeldbo, and V. R. Eshleman, "Models for the Atmo- 
sphere of Mars, Based on the Mariner IV Occultation Experiment," 

J.G.R. , 71, No. 9, 2307-2316, 1966. 

8. Evans, J. V. and G. H. Pettengill, Memorandum to JPL Advisory Group 

on Radio Experiments in Space; Re: Bistatic-Radar Studies of the 

Moon, 23 January, 1961. 

9 Eshleman, V. R. et al., "Bistatic-Radar Astronomy," Chapter 6 of Summer 
Study 1962 Review for Space Science Board NAS/NRC-1079, Iowa City, 

Iowa . 

10. Private communication by V. R. Eshleman 

June 1960: Meeting on Electromagnetic Studies of the Moon and 

Planets Using Space Probes," at Jet Propulsion Laboratory, Pasadena, 
California. 

June 24, 1960: "Conference on Planetary Atmosphere," Space Sciences 

Board of National Academy of Sciences, Arcadia, California. 

December 15, 1960 and February 2, 1961: "Conference on Planetary 

Atmospheres," Space Sciences Board of National Academy of Sciences, 
California Institute of Technology, Pasadena, California. 

February, 1961: Meeting of Ionospheric Physics Subcommittee of 

Space Sciences Steering Committee, NASA, Stanford, California. 

11. Tyler, G. L. et al.: "Bistatic-Radar Detection of Scattering Center 
on the Surface of the Moon with Lunar Orbiter," In press, April 1967. 

12. Helstrom, C. W., Statistical Theory of Signal Detection, Pergamon 
Press, New York, 1960. 

13. Ibid , cf. Chapter II, Section 5. 

14. Jackson, J. D., Classical Electrodynamics , Chapter 6, John Wiley and 
Sons, New York, 1962. 


107 


SEL -67-042 


15 . 


Born, M. and E. Wolf, Principles of Optics , Chapter 10, MacMillan, 

New York, 1964. 

16. Ibid , Chapter 4. 

17. Brouwer, D. and G. M. Clemence, Methods of Celestial Mechanics, 

Academic Press, New York, 1961. 

18. Beckman, P. and A. Spizzichino, The Scattering of Electromagnetic 
Waves from Rough Surfaces , Chapter 7, Pergamon and MacMillan, New 
York, 1963. 

19. Turin, G. L., M An Introduction to Matched Filters,” IRE Trans, on 
Information Theory, 311-329, June 1960. 

20. Sherwin, C. W. et al . , "Some Early Developments in Synthetic 
Aperture Radar Systems," IRE Trans. ME, MIL-6, 111-115, April, 1962. 

21. McCord, H. L., "The Equivalence among Three Approaches To Deriving 
Synthetic Aperture Array Patterns and Analyzing Processing Tech- 
niques," IRE Trans. ME, MIL-6, 116-119, April, 1962. 

22. Cutrona, L. J., "A Comparison of Techniques for Achieving Fine 
Azimuth Resolution," IRE Trans. ME, MIL-6, 119-121, April, 1962. 

23. Heimiller, R. C., "Theory and Evaluation of Gain Patterns of 
Synthetic Arrays," IRE Trans. ME, MIL-6, 122-129, April, 1962. 

24. Tyler, G. L. , "The Bistatic, Continuous-Wave Radar Method for the 
Study of Planetary Surfaces," J. Geophys . Res., 71, 6, 1559-1568, 

March 15, 1966. 

25. Gabor, D., "A New Microscopy Principle," Nature, 161, 771, 1949. 

26. Leith, E. N. and J. Upatnicks, "Reconstructed Wavefronts and 
Communication Theory," J. Optical Soc. Am., 52, 10, 1123-1130, 

October, 1962 . 

27. Courant, R. and D. Hilbert, Methods of Mathematical Physics , Vol. I, 

Interscience Press, New York, 1965, p. 138. 

28. Helstrom, C., Statistical Theory of Signal Detection , cf . Chap- 
ter IV. C., Pergamon Press, New York, 1960. 

29. Hofstetter, E. M. , "Some Results on the Stochastic Signal Parameter 
Estimation Problem," IEEE Trans. Information Theory, IT-11, 422-429, 
July, 1965. 

30. Helstrom, C., Statistical Theory of Signal Detection , cf . Chapter XI, 
Pergamon Press, New York, 1960. 

31. Middleton, D., Introduction to Statistical Communication Theory, 
McGraw-Hill, New York, 1960. 

32. Price, R. and P. E. Green, "Signal Processing in Radar Astronomy 
Communication via Fluctuating Multipath Media," MIT Lincoln Laboratory 
TR No. 234, 6 Oct. 1960. 

33. Middleton, D., "On New Classes of Matched Filters and Generalizations 
of the Matched Filter Concept," IRE Trans, on Information Theory, 

IT-6, 349-360, June 1960. 


SEL-67-042 


108 



34. Yakovin, A. A., "Motion, Rotation, and Figure of the Moon," Ch, 1, 

From The Moon: A Russian View, ed. A. V. Markov, Univ. of Chicago 

Press, Chicago, Illinois, 1960. 

35. Tyler, G. L. et al . , "Bistatic-Radar Detection of Scattering Centers 
on the Moon Using Lunar Orbiter I," in press, March 1967. 


109 


SEL-67-042 


Security Classification 


DOCUMENT CONTROL DATA . R & D 


(Security elatsltlemlton ot title, botly of •bslrmcl and Indexing annotation muni be ./irered when the q W „h , aDO „ u cU* SIUed) 
ORIGINATING ACTIVITY fCoipor • •Uthot) l 2 .. REPORT SECURITY CLASSIFICATION 

Stanford Electronics Laboratories Unclassified 

Stanford University, Stanford, California 1 2b GROUP ~ ’ *“ 


3. REPORT TITLE 


BISTATIC-RADAR IMAGING AND MEASUREMENT TECHNIQUES FOR THE STUDY OF PLANETARY SURFACES 


4. DESCRIPTIVE NOTES (Type oi report and inclusive dates) 

Scientific Report No. 19 


5. au THOR(S) (First name, middle initial, last name) — ■ . ... 

George Leonard Tyler, Jr. 


6- REPORT DATE 7a. TOTAL NO. OF PAGES ( 76 . NO. OF REFS 


Ba. CONTRACT OR GRANT NO. 


Grant NsG-377 


b. PROJECT NO. 


9fl. ORIGINATOR'S REPORT NUMBER(S) 


Scientific Report No. 19 


9b. OTHER REPORT NOIS) (Any other numbers that nuty be assigned 
this report) 

SEL-67-042 




13. ABSTRACT 


12- SPONSORING MILITARY ACTIVITY 


National Aeronautics and Space 

Administration 


The problem of imaging a planetary surface is considered from the point of view of 
continuous wave bistatic-radar in which transmissions originating on the Earth (or on 
a spacecraft) would be received on a spacecraft (or on the Earth) after reflection froir 
a planetary surface. 

An electromagnetic wave reflected by a planet may be considered as a superposition 
of the waves from elementary scatterers on the surface. Using a superposition model 
the ordinary radar brightness distribution, the differential radar cross-section, the 
polar scattering diagram, and the radar albedo of a surface are all related measures 
of the local surface properties. 

The surface brightness distribution may be determined by cross-correlating the 
scattered fields (as measured by a spacecraft along some fraction of its trajectory) 
with the expected signal from each point on the surface. The predicted azimuthal res- 
olution in wavelengths is inversely proportional to the angle subtended at the target 
point by the fraction of the trajectory over which the data are taken; in range the 
resolution in wavelengths is inversely proportional to the square of the same angle. 

The feasibility of the method depends on the use of the illuminating wave on board the 
spacecraft as a frequency reference to achieve the requisite stability. 

Physical analogs of the process exist as modifications of holograms and synthetic 
antenna arrays. An additional analog is that of a bank of tracking filters, each of 
which is adjusted to receive the signals from a separate portion of the surface. 

The maximum-likelihood estimator for the brightness of a specific scattering area, 
in the presence of white Gaussian noise at low input signal-to-noise ratios, is a 


1473 


UNCLASSIFIED 

Security Classification 








UNCLASSIFIED 
Security Classification 


1 4. 

K EV WORDS 

RADAR 


BI STATIC -RADAR 


IMAGING 


PLANETARY SURFACES 



Abstract (cont 

Hilbert quadratic form with the expected auto- 
correlation function of the signal from the 
scattering area as a kernel. For maximum res- 
olution (minimum area) this is equivalent to the 
cross-correlation method for obtaining the bright- 
ness distribution. Furthermore, this estimate is 
efficient in the sense of being unbiased and of 
having minimum variance, for the low signal-to- 
noise ratio case. 

The estimator may be realized with well-known 
forms of time-varying or time-invariant filters, 
or with correlators. A new realization is that of 
a hologram scanned with a properly weighted illu- 
minating wave . 

It is concluded that bistatic-radars operatin 
between a ground station and a vehicle in space 
constitute a powerful technique with which to map 
and to study planetary surfaces. 



UNCLASS I FI 


Security Claaeification 










